A SIMULATION MODEL FOE STUDYING EFFECTS OF POLLUTION AND 
FRESHWATER INFLOW ON SECONDARY PRODUCTIVITY IN AN ECOSYSTEM 



fey 


N75- 1 3405^ 

Uncl'as ' v 
G3/45 053 1.9' 

A thesis submitted to the Graduate Faculty of 
North Carolina State University at Raleigh 
• in partial fulfillment of the 
requirements for the Degree of 
’• Doctor of Philosophy 

DEPARTMENT OF MARINE- SCIENCES 


RALEIGH 

1 9 T k 


ROBERT WARD JOHNSON 


; (N ASA»TM-X-72169) A SIMULATION MODEL FOR 
{STUDYING EFFECTS OF POLLUTION AND 
(FRESHWATER INFLOW ON SECONDARY 
] PRODUCTIVITY IN AN ECOSYSTEM Ph.D. 

'Thesis - (NASA) - *CSCL 13B 




L 7 


APPROVED BY: 


Co-Chairman of Advisory Committee 




- Reproduced by 


l - NATIONAL TECHNICAL 
• INFORMATION SERVICE 

U S Department of Cemmaree 
__ „ .Springlieid, VA. 22151 _ 







ABSTRACT 


JOHNSON, ROBERT WARD. A Simulation Model for Studying Effects' .of 
Pollution and Freshwater Inflow on Secondary Productivity in an 
Ecosystem. 

Operations research methodology is used to develop a mathematical 
model of the Galveston Bay, Texas, ecosystem. Secondary productivity, . 
measured by harvestable species (such as fish, crabs, and shrimp), is 
evaluated in terms of man-related and controllable factors, such as 
quantity and quality of inlet fresh’ water and pollutants. The simula- 
tion type model uses information from an existing physical parameters 
model as well as pertinent biological measurements.. 

One of the purposes of the model is to provide predictive informa- 
tion of value to those responsible for estuarine management. Results 
are of major benefit in pollution control and fisheries management in 
estuarine systems, particularly those which have migrating species, 
which include fish (menhaden, trout, bass, croaker for example), shrimp, 
and crabs. Another objective of the study is to identify those biologi- 
cal, chemical, and physical parameters that should be measured in order 
to develop models for similar ecosystems. 

The Galveston Bay, Texas, is a highly productive temperate-zone 
ecosystem that has been subjected to man-related stresses. There are 
extensive hisotrical biological data and analyses available as well as 
chemical and physical parameters models of the ecosystem. This in- 
depth information is highly desirable for the modeling of an ecosystem. 
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INTRODUCTION 


Estuarine systems , due to their unique locations at the interface 
of rivers and seas, are under heavy user . and development pressures. - 
These highly competitive forces include .economic , political , social, . 
and recreational potentials and place high priority on certain estuarine 
characteristics that are important for man-related uses for commercial 
(manufacturing and shipping), residential, and recreational purposes, 
as well as "biological productivity. Concurrent incompatible uses in • 
close proximity lead to "pollution" conditions; that is, a user benefit 
suffering due to another activity. 

Generally, the concept of pollution involves a man-related activity 
that disturbs the ''natural" system; or more specifically, the discharge 
of nutrients or other energy (heated water) and/or toxic materials. 

Other forms of pollution (or something that perturbates the system) 
involve flow characteristics of the system, consequently affecting the 
environmental conditions such as temperature and salinity. These may 
be due to dredging in the estuary, or upstream river flow changes such 
as damming and impoundments. These factors affect the estuarine system, 
whether or not combined with nutrient or toxic discharges (Gunter 1950, 
196l; Gunter, Christmas and Kellibrew 1964; Copeland 1966; Cooper 1970; 
Copeland and Bechtel 1971). 

Each of the users which may potentially benefit from the character- 
istics of an estuarine system should have an input and consequently 
some priority to obtain estuarine benefits. In any event, use of 
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estuarine systems is the result of management, whether formal or 
informal. As competive uses grow it will be necessary to apply a 
more formalized management • system, which requires analyses that will . 
provide information so the best decisions may be made. This general 
approach is within the scope of the systems analysis ‘(Van Dyne 19 ^ 9 ; ■ 
Dale 1970 ; E.-P. Odum 1971;. Patten 1971), and requires that quantita- 
tive information be available on which to base the analysis and subse- 
quent inputs to. management , who make the decisions. 

One of the areas currently notable for the lack of quantitative 
information to use as inputs to the decision-making process is associ- 
ated with the biological complex of an estuary and the resultant 

6 

changes in the system as various exogenous changes are made, whether 
natural or man-caused. 

Notable and significant groundwork, some theoretical , some .para- 
metric, has been initiated by Patten (1959) » H. T; Odum ( 1967) and . 
Paulik (1971). Computerized models have been developed for 'certain 
specific cases, such as the salmon industry, by Paulik (1967)* 

On an ecosystem conceptual basis. Watt (1966, 1968), Van Dyne 
(1969), Koval (1971) and Patten (1971) have developed field, insect, 
and aquatic ecosystem models to demonstrate systems analysis as an 
applicable quantitative tool. Williams (1971), in his computer simu- 
lation of Linderman’s (19^2) data from Cedar Bog Lake, indicates the 
lack of data available for current modeling procedures. 

One of the areas probably receiving too little attention in past 
modeling efforts have been the necessary steps of following through on 
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the modeling process by making predictions (using an initial model) , 
gathering current data, comparing the predicted to- observed data, 
making any necessary model corrections, then repeating the process until 
suitable agreement is obtained (Watt 1968). Since a system reaction to 
exogenous changes should be the same, historical or new data may be used 
in the above modeling processes. , 

Fisheries studies in estuarine systems have concentrated on: 1) •. 

migration (Gunter 1950; Copeland 1965; Copeland and Fruh 1970), 2) feeding 
and food availability (Darnell 1958, '1959* 196l; Hea-ld '1971; W. E. , 

Odum 1971), and 3) pollution effects on respiration and growth (Wohlschlag 
1972). Results of these studies provide qualitative and short term 
cause and effect analyses within an estuary. In addition, the concept 
of diversity index (Copeland and Fruh 1970i Copeland and Bechtel 1971) > 
and its relation to pollution have aided in the quantitative assessment 
of pollution on the biological communities. 

International and open sea fisheries studies have been concerned 

. , . * a 

with growth and population dynamics of the fish of commercial interest. 
Pertinent information on growth characteristics are developed by von 
Bertalanffy (1938), Parker and Larkin (1959)* Ivlev (1966)., Ursin (1967), 
and Nickolski (1969). Most of these concepts are based on theoretical 
as well as sampling and commercial catch studies, and consequently 
have broad application. 

Environmental stresses are commonplace in natural systems and are 
a large factor in determining communities and populations (Slobodkin 
i960, 1962, 1967; Margalef 1963, 1968; Odum, Copeland, and McMahon 1969)* 
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Estuarine systems, particularly those in the temperate regions , have 

■ ' .1 ' 

seasonal stress patterns that dominate the natural biological communi- 
ties through control of these environmental stress factors. The 
seasonal peak of energy flow down the rivers in the spring leads to 
the well-known migrations of commercially important fish, shrimp, and 
crabs to the estuaries where not only is there bountious food, but 
essentially predator-free nursery areas. Rapid growths occur, followed 
by outward migrations due to organism physiological preference and envi- ' 
ronmental changes of temperature and salinity (Copeland 1965; Copeland 
and Truitt 1966). 

In recent years , as industrial and population growths have increased 
along rivers and estuarine areas, a new set of stresses have been levied 
on the biological populations. An in-depth study of the effects of 
.thermal, nutrient, and toxic effluents -was . reported by Brett (1957) on 
the rivers of British Columbia. The basic approach was to evaluate the 
effluents as causes of stress (with stress being defined as' a state 
under which chances for survival are reduced) taking into account the 
entire life span of the organism as well as short-term effects. Cronin 
and Flemer (1967) also evaluated the effects of pollution on energy 
transfer in coastal environments and include chemical growth inhibitors 
as a significant factor. 

Toxic pollution causes stress conditions that are indiscriminate 
in character; i.e., affects each member of the group (as opposed to 
discriminate stress which affects individuals single, but not the group 
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as a whole, like predation, individual, parasitism, trapping, etc.) and 
may be lethal, limiting, inhibiting, or loading (Brett 1957)* These 
three latter stress conditions and their effects on growth, survival, 
metabolism, and population dynamics have also been 'investigated by 
Steed and Copeland (1967), Wohlschlag and Cameron (1967), Mount .(1968), 
Wohlschlag, Cameron and Cech (1968), Copeland and Wohlschlag (1971)* 
Kloth and Wohlschlag (1972), and Wohlschlag (1972). There are many 
factors associated with sublethal stress conditions that are still not 
clearly understood, but it is evident that the decreased metabolism 
due to toxic and temperature effects (Warren and Davis 1966 ; Copeland 
and Wohlschlag 1971; Wohlschlag 1972) leads to decreased growth rates 
and subsequent effects on the ability of species to compete and survive 
It is the purpose of this effort to develop a biologically sound 
computerized simulation model of the biological energy flow through an 
estuarine system, specifically the Galveston Bay, Texas. Fish, shrimp, 
suid other organism growth characteristics will be based On logistic 
growth patterns as recommended by von Bertalanffy (1938), Ursin (1967), 
and Patten (1971)* with limiting conditions imposed by food supply, 
migrations,' and stresses (toxic or environmental). Exogenous inputs 
will include "natural" variations such ds seasonal immigration and , 
emigration, salinity, and food supply patterns. In addition , man- 
controllable (through management), factors related to river freshwater 
flow manipulation and pollution effluent will be considered. System 
outputs will be biomass levels for the organisms (fish, shrimp, etc..) 
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in the bay and emigration, which will be taken as a measure of the 
productivity of the bay. 

Validation of the model will be approached by 1) comparing organism 
food consumption to that available in the estuary and 2) comparing pre- 
dicted results from the model to results obtained from analogous field 
studies. In the latter case the model is used to investigate the effects 
on fisheries, productivity due to changes of pollution and/or further 
decrease in freshwater (e.g. , due to damming of the Trinity River). 

A further objective of this model and analysis effort will be to 
identify areas of future research and/or data needed for effective and 
efficient pollution control and management of estuarine systems. 
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THE GALVESTON. BAY 

Physical , Chemical , and Biological Characteristics 

Physical, chemical, and biological characteristics of the Galveston 
Bay, Texas, estuary have been studied by Cooper (1970), Copeland and 
Fruh (1970), and Armstrong and Hinson (1973) • Major effects in the 
system are due to changes in freshwater inflow and pollution inputs . 
These affect the. primary characteristics of salinity an.d total nitro- 
gen (which is used as a measure of pollution load in this study). Aver- 
age annual salinity and total nitrogen distributions for 19&9 (Copeland 
and Fruh, 1970) are. shown in figures 1 and 2, respectively. 

Copeland and Fruh (1970) reported percent sources of water 'at each' 
of their sampling stations in Galveston Bay, obtained from a low- flow 
conservative model. Source water fractions and yearly • average salinity 
and total nitrogen values are listed in table I. Salinity and total 
nitrogen values for subsequent years were determined from source water 
changes on a station by station basis. Estuary levels were determined 
by arithmetically averaging the station values. 

Armstrong and Hinson (1973) investigated freshwater inflow quan- 
tities and waste discharges from the major tributaries of the Galveston 
Bay. These. values were grouped into the same sources as used by Cope- 
land arid Fruh (1970) for each of the stations for subsequent analysis of 

t 

changes. Ten year water and three year (1969-71) "waste average discharge 
rates were 




Figure 1.- Isohalines for the Galveston Bay system plotted from mean 

annual salinity (ppt) at each station. Data from collection 
cruises during February., April, July and October 1969 (from 
Copeland and Fruh , 197 0 ) . 
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Figure 2.- Mean annual concentration gradients of total nitrogen (mg/l) 
for the Galveston Bay system. Data taken from monthly values 
of the 1969 Bay Sampling Program (from Copeland and Fruh, 

197 0 ) . 
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Table I: Galveston Bay station water sources and yearly average salinity (ppt) 
and total nitrogen (mg /i) from Copeland and Fruh, 1970. 




Source Water 




Station 

Bol. Roads 

T. River 

HSC 

Other 

Sal, ppt 

T NIT, mg/1 

t.oo 

.74 

07 

~~.lV’ 

.05' '• 

27.00 ' 

0.00 

2.00 

.50 

0.00 

.10 

.40 

25.22 

.75 

3.00 

.50 

0.00 

.50 

0.00 

. 22.80 

■ .60 

4.00 

.40 

.10 

.50 

0.00 

17.00 

1.00 

5.00 

.25 

.10 

.65 

0.00 

15.50 

1.20 

12.00 

.70 

0.30 

0.00 

.30 

22.50 

. 80 

13.00 

.80 

0.00 

0.00 

.20 

29.00 

.40 

14.00 

. 75 

0.00 

0.00 

.25 

24.00 - 

.60 

15.00 

.47 

.16 

.33 

.04 

19.00 

.92 

16. 00 

.49 

.09 

.18 

.24 

21.30 

.82 

17.00 

.50 

.09 

.18 

.23 

19.50 

U05 

lfl.00 

.43 

.16 

.37 

.04 

13.50 

.89 

l 9. 00 

.29 

.16 

.55 

.01 

1 1.50 

1.15 

20.00 

.24 

. 14 

.61 

.05 

11.50 

1.70 

21.00 

.13 

.09 

.78 

-0.00 

14.00 

3.00 

22.00 

.26 

.15 

.58 

.01 

14.80 

. 1.91 

23.00 

.24 

.42 

.36 

0.00 

12.80 

1.90 

24.00. 

.20 

.54 

.28 

0.00 " 

, '• 4.50'. 

1.00 

25.00 

.16 

• 63 

.22 

0.00 

" - 2.50 

l . 00 

26.00 

.21 

.52 

.29 

• 0.00 

7.40 

.70 

27.00 

2 r.oo 

.23 

.48 

•31 

0.00 

7.00 

.90 

.31 

.34 

-AT.. 

0.00 

9.00. 

.95 

29.00 

.37 

.27 

.37 

0.00 

11.00 

.90 

30.00 

.35 

.27 

.38 

0.00 

11.00 

.90 

3L.00 

.5 1 

.09 - 

.17 

.23 

2 2.50 

.70 

32.00 

.49 

.09 

.17 

.25 

21.00 

.70 

36.00 

.14 

0.00 

• 86 

0.00 

13.00 

5.00 
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Source 

Trinity Houston 
River Ship Channel Other Total 

Water Discharge, 1000 CFS 5-58 10. 2.79 18-90 

Waste Discharge, million 29*90 153.2 52.20 235-3 

pounds of BOD/yr. 

Analysis of rainfall data (U.S. Weather Bureau) indicates that 1969 
and the preceding two years were ’’average' 1 . 

From these water and waste discharge information and the seasonal 
detritus curves of Heald (1971) for. a South Florida Bay a detritus input 
curve was determined for the Galveston Bay. Significant considerations 
were that l) biological oxygen demand (BOD) reported by Armstrong and 
Hinson (1973) inherently included waste products as well as detritus 
(organic particulates); 2) due to similar ecosystems (having many of 
the same species) marsh and submerged grasses would tend to have similar 
seasonal cycles, even if displaced by several months; and 3) maximum 
and minimum value ratios of detritus densities would be approximately 
the same for the two ecosystems. The resultant detritus (organic par- 
ticulates) input curve to the Galveston Bay for the calendar year 1969 
is shown in figure 3. Detritus loads for subsequent years are discussed 


in Appendix A. 
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Results of biological' 'sampling in the Galveston Bay in 1969 (Copeland 
and Fruh 1970 ) are summarized in figure 4 . Paired curves are shown for the 
zooplankton counts and other consumer groups biomass (consumer groups will 

• •'I 

be defined in the following section). Biomass values are the sums from 
the individual stations, and will be reference levels in this study, 
rather than estimating total ecosystem biomass values. 

The Galveston Bay Ecosystem 

Temperate zone ecosystem characteristics are dominated by an annual 
seasonal cycle that is, in general, controlled by weather (Chin, 1961). 
Processes in the ecosystem are related to energy sources (food) and 
migrating consumers with physiological adaptations that allow them to 
effectively compete for the available foods. 

■ In the Galveston Bay ecosystem freshwater flow from the feeding rivers 
brings in large quantities of organic particulates (detritus) and dis- 
solved nutrients which serve as energy sources for the base of the feed- 
ing chain. Dissolved nutrients are necessary for growth of, marsh grasses 
( Spartina spp ), fixed bottom plants (turtle grass, Thalassia testudinum , 
for example) and small floating plants (phytoplankton). These materials 
are grazed or filtered from the water by small animals such as zooplank- 
ton, herbivores (shad, Dor o soma cepedianum and menhaden, Brevoortia par- 
tronus ) and omnivores (shrimp, Penaeus spp and crabs, Callinectes sapidus ) « 
These small animals are in turn consumed by larger carnivore species 
(Atlantic croaker, Micropogon undulatus . Anchovy, Anchoa ini t chilli and 
Trout, Cynoscion arenarius ). A generalized energy flow diagram for the 
Galveston Bay is shown in figure 5.' 



Log biomass, grams (log count for Consumer Group 1) 



Energy 

sources Foods 


Foods and 

consumers Consumers 


Solar 


Dissolved 

nutrients 


Detritus 


Vascular 
plant mtl. 


Zooplankton 
Consumer Gp. 1 


Herbivores 
Consumer Gp. 2 


Mixing 

mi 


2S « 

Omnivores 



Consumer Gp. 3 


V 



Prim carnivores 
Consumer Gp. 4 


Mid. carnivores 
Consumer Gp. 5 



Top carnivores 
Consumer Gp. 6 


Figure 5.- Biological- energy flow in Galveston Bay ecosystem. 
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One ecosystem characteristic is that certain functions are performed 
"by one or more species, either simultaneously or over a period of time. 
For this reason, the biological species in the Galveston Bay may be 
grouped based on similarity of function and feeding characteristics'. In 
this study consumer groups have been organized based on consuming habits 


food preferences.. 

Consumer groups and typical 

species in them are: 

Consumer 

Group 

No. 

Consumer Group 

Typical Members 

1 

Zooplankton 


2 

Herbivores 

Menhaden 

3 

Omnivores 

Shrimp 

h 

Primary Carnivores 

Atlantic Croaker 

5 

Middle Carnivores 

Anchovy 

6 . :• 

Top Carnivores 

Trout 


The above groupings are based on dominant characteristics in the 
first year or period of maximum rate of growth in the ecosystem. Adults 
may not consume the same foods as the young of the same .specie s* however, 
this is not a limitation since consumer groups and shifting of consuming 
habits are included in the model. 

Feeding habits of estuarine species have been studied extensively by 
Darnell (1958, 1961) and W. E. Odum (1971). Based on their results foods 
curves have been developed for the six consumer groups, figures 6-11. 

Productivity (i.e., growth of biomass in the estuary followed by 
harvesting or catches) either in the estuary or after it has left the 
estuary, is one of the uses of an estuarine system. It provides food, 



Fraction energy 






Fraction energy 
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employment and other economic benefits and is one measure of the value 
of an estuary. 

Productivity in the estuarine system as indicated previously is 
determined by immigration, growth, and emigration. Thus, the estuary 
meets a specific need for the species of interest during some period of 
their lifetime. Immigration and emigration are largely natural phenom- 
ena that represent an adaptation of particular species to' the total en- 
vironment in which they live and are not directly controllable or manag- 
able by man. (Perhaps this is also due to the fact that the total system, 
which includes the oceans, is too big!) On the other hand, the growth 
phase of estuarine organisms is highly affected by man's activities, 
particularly pollution due to waste discharges and manipulation of water 
flows into of within the estuarine system (Odum, Copeland, and McMahon 

1969 ). . 

Factors identified as primary in the growth of estuarine species 

are: 

1) Food and consumer densities; 

2) Environmental effects, including pollution; and 

3) Distribution effects, which are related to geological parameters. 
Other parameters, such as temperature, are obviously important, (£.£. , 
Gunter, 1957 ) but do not appear to be controlling factors in changes in 
the Galveston Bay at this time. In any event, it appears that their 
relation to growth per se are secondary compared to the three factors 
listed above. 

Food and consumer densities and their effect on consumer growth 
rates have been studied by Brocksen, Davis and Warren (1970). In the 
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relatively confined basins and lakes in which their studies were con- 
ducted on sockeye salmon ( Oncorhynchus nerka ) there was a linear rela-. 
tion between consumer biomass and growth rate, figure 12. Extrapolation 
to a zero growth rate at some consumer density is questionable, and, as 
they pointed out, even in the lowest food producing areas, they did not 
measure a zero growth rate. Consumer growth rate as a function of food 
density was also a part of the above study. In this case food preference 
of the salmon was very specific to zooplankton and salmon growth rates 
had a direct correlation to zooplankton densities over a wide range, 
figure 13. 

Effects of pollution, whether nutrients or toxic materials, has been 
well established and has significant, effects not only on specific species 
but on communities as well (Copeland, 1966; Steed and Copeland, 1967; and 
Wohlschlag and Cameron, 1972). However, quantitative relationships over 
the period of an organisms or species life cycle has not. been established. 
Alderdice and Brett (1957) and Brett (1957) investigated the effects of 
kraft paper-mill wastes and hydroelectric power plants on the growth and 
survival potential for salmon during their migrations. One of their sig- 
nificant. concepts was that relatively small increases in stress (due to 
pollution or environmental factors) could lead to significant loss of 
competitiveness and possible elimination. Steed and Copeland (1967) 
reached similar conclusions from studies on pinfish ( Lagodon rhomboides ) 
using petroleum waste effluents. Wohlschlag (1972) studied the effects 
on metabolism of the relatively polluted Galveston Bay waters compared 
to those in Aransas Bay (collected on incoming tide thus essentially 



Growth rate (K x 10 ) 


— A— Production 
— O — Growth rate 



Figure 12..- Consumer growth rate as a function of consumer biomass 
(taken from figure 2, Brockson, Davis, and Warren, 1970) 
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Figure 13.- Consumer growth rate as a function of food availability 
(taken from Brocksen, Davis and Warren, 1970).. 
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pollution free) using striped mullet ( Mug 11 cephalus ) . His results 
indicated about a 10 % reduction in metabolism and probably at least 
that "much reduction in growth. 

Freshwater inflow into Texas Coastal Bays is a significant factor 
in their productivity. In addition, as observed by Hildebrand and 
Gunter (1953), Copeland (1966), and Armstrong and Hinson (1973), the 
effects are species specific and, on commercial catches (such as shrimp) 
may have a one to two year lag. Copeland (1966) studied the effects of 
freshwater inflow as a function of bay or estuarine size and location; 
Armstrong and Hinson (1973) subsequently refered to this as displacement 
rate. In general in .a specific estuary, such as the Galveston Bay, 
total commercial catch increased with decreased freshwater input (at 
least to the point of a significant ecosystem shift); however, the 
decreased freshwater flow led to higher production of f infish species 
such as the Atlantic Croaker at the expense of the more economically 
desirable shrimp and crabs. Armstrong indicated this may be due to 
decreased spawning and feeding areas as a result of higher salinities 
in the bordering marsh areas. Data from the above studies are discussed 
in section III on model calibration. 

It is interesting that, largely due to the magnitude of the prob- 
lem, it has been only recently that studies of multiple effects and 
their interactions have been initiated and analyzed (e.g., Alderdice 
1963, 1972). However, we would intuitively recognize that an environ- 
mental effect is not constant but varies as a result of other parameters 
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This condition has "been well established for salinity-pollution (Cope- 
land and Fruh, 1970) and for temperature-growth (Kinne 1965? 1967); 
however the effect has not been quantitatively established in a natural 
growth situation, as is being evaluated in this study. Identified 
growth effects used in this study have been reported from prior studies 
of them as independent 'factors . Interactions may be included in some 
of the effects reported; however, they are probably of secondary 
importance, as for temperature. 
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SIMULATION MODEL OF THE GALVESTON BAY 

A continuous simulation model format (Forrester 1961 ) is used for 
the Galveston Bay ecosystem. Independent variables are exogenous changes 
in freshwater and waste discharge to the ecosystem. Dependent variables, 
or outputs, are biomass levels of six identified consumer groups. Analy- 
tical and empirical relations are used to define and relate physical, 
chemical, and biological characteristics of the ecosystem. 

This section provides a functional description of the model. A 
detailed description, , including equations, is in Appendix A and a 
program print-out in Appendix B. For clarification, it may be useful 
to review the previous section describing the ecosystem. 

The model 

In an estuarine ecosystem, as discussed previously, the dominant 
cycle is the seasonal calendar year. In the model, the calendar year 
is divided into 26 two week periods (designated as I = 1,,,.,26). These 
periods were short enough that rate changes within the period are 
insignificant. 

Each of the six consumer groups (designated as L - 1,...,6) is 

* 

phased into the yearly cycle, but their own cycle is different from 
the others, as shown in consumer biomass curves developed from sampled and 
historical data for the initial model (baseline) year, figure i*. Each con 
sumer group's seasonal cycle (periods designated by M) starts with.M = 1 
defined as the period when a consumer group’s biomass is at a minimum 
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in the ecosystem (e.g. , M = 1 in period I = 21 for consumer group 4). 

The biomass in the estuary at M - 1 is assumed to be residual, and is 
therefore designated as ."adult” (as contrasted to this year's young or 
the immigration into the ecosystem). The seasonal cycle and its rela- 
tionship to the calendar year is shown schematically in curve "e" of., 
figure 14.- . • 

Immigration starts in period M = 1 and continues over 3-A model 
periods (6-8 weeks). -Functionally, it is taken to be a sinusoidal 
shaped curve (positive l80°) and is shown schematically as curve ”a" 
of figure lU„ Immigration is a constant from year to year, independent 
of ecosystem variations (Caillouet and Baxter 1973; Copeland 1973). 

Immigrated (larval and post-larval) organisms have very high growth 
rates decreasing with increasing organism size (Paloheimo and Dickie 
1965; Patten 1971). A decreasing value exponential function is used 
for the base-line year to describe the growth rate over a consumer 
group's year, decreasing to a value of 0 for adults. (e,.g,. , beginning 
of the next year). Base-line year growth rates are shown; schematically 
as curves "b" and "c” for immigrated and adult organisms, respectively. 
Note that in the base-line year all growth rate effects due to changes 
in exogenous variables have, by definition, values of 0 (or a multi- 
plier of 1.0). In subsequent years, year to year changes in the exoge- 
nous variables are defined in terms of growth rate change ratios , which 
are used to determine new values for net growth rates. 

After high growth rates in the estuary, consumer organisms emigrate 
from the ecosystem - in this case, primarily to the Gulf of Mexico. 
Emigration is shown schematically as curve ”d” in figure l4v 
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Figure 14.- Development of baseline year biomass , curve for consumer 
group 4 . 
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Consumer group h biomass for the base-line year is shown schema- 
tically as curve "e" in figure 1^. Biomass curves (counts for zooplank- 
ton) for the other five .consumer groups are independently developed .in 
an analogous manner.- These six curves are correlated with the calendar 
year periods (i), which are used for period- identification after the 
base line year. 

In the model, iterative calculations are made period by period. A 
period calculation consists of l) determining the net period growth 
rate, which is the product of the prior year net growth rate and current 
year change ratios due to consumer and food densities and exogenous 
variables (in the baseline year the net growth rate is taken from the 
exponential curve); 2) to the biomass at the beginning of the period 
adding immigration and subtracting emigration; and 3) the total is 
multiplied by the net growth rate to obtain period biomass which is-, 
also the biomass at the beginning of the next period. 

Model stability is aided by built-in safeguards which act as 
negative feedback; first, if environmental conditions remain the' same, 

biomass curves will repeat the previous year's values, except for time- 

<; ’ 

lag effects; second, in the food and consumer density function the 
growth rate ratio is an inverse function of consumer density compared 
to the baseline year value (in effect this implies an upper limit on 
consumer biomass in the ecosystem) (equation 8, Appendix A); and third, 
model calibrations and examples are based on wide ranges of freshwater 
and pollution inputs that occurred over about a 20 year period (sections 
II, l,b and IV). . 
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The Baseline Year 

In 1969 a comprehensive sampling program was accomplished as part 
of the Galveston Bay Program (Copeland and Fruh, 1970 ). Results of 
this sampling program have been used to determine numerical parameters 
for seasonal changes in consumer group biomass levels for the Galveston 
Bay pollution effects model (i.e., biomass curves expressed in model 
language). For this base-line year, all growth change ratios are set 
equal to 1 . 0 . This allows the development of growth rates for the 
conditions that existed during the base-line year. 

As developed in the previous section, biomass curves for each of 
the consumer groups includes immigration, growth in the ecosystem, 
and emigration. Typical of lower temperate region ecosystems, year- 
round populations remain in the estuary. In the model- the biomass in, 
the ecosystem at the minimum biomass of a consumer group is taken as 
"adult" or the year-round population. As noted previously in the model' 
section, this period of minimum biomass is the start of that consumer 
group's seasonal year (M = 1 for that consumer group). Pertinent 
events in the development of the base-line year biomass curve and numer- 
ical values for consumer group h are (see figure l4): 


Period 

Event* ^1 M.** 


Start of Immigration 

21 

1 

End of Immigration 

25 

' 5 

Start of Emigration 

. h 

10 

End of Emigration 

21 

1 


#See Appendix A. for model detail. 
#*M relates to baseline year only 
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Numerical Values* 


Constant of Immigration 

1600 

Factor of Emigration 

.20 

Growth Rate Factor 

1.20 


#See Appendix A for model detail. 

The model biomass curve is shown in figure 15 . Values from the 
sampling program (Copeland and Fruh 1970) > are also shown. 


Model Calibration 

Model calibration, or' development of numerical equations to describe 
changes in. growth rates as a result of changes in food and consumer 
densities and exogenous variables, is based on a number of independent 
investigations, - each of limited scope. In the model, as in the studies, 
factors other than those being evaluated will.be held constant during 
that phase of model calibration. Effects on consumer group growth rates 
will be evaluated in the following order: 

1) Food and consumer densities; 

2) Environmental effects; and 

3) Distribution effects. 

Growth rate effects due to consumer and food densities were studied 
by Brocksen, Warren and Davis (1970) as discussed previously, figures 
12 and 13* Specific model equations are developed in the Appendix. 




Calendar year periods (I) 


Figure 15.- Resulting year model biomass, curve and sample data for. consumer group 1 
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Environmental effects on growth - primarily due to pollution 
effects - have been investigated by Wohlschlag (1972) as discussed 
previously. Growth effects are assumed to vary linearly with the pol- 
lution parameter. 

Distribution effects in the Galveston Bay are due primarily to 
changes in freshwater inflow. Investigations by Copeland (1966) and 
Armstrong and Hinson (1973) have evaluated shrimp and total biomass 
productivity as influenced by freshwater inflow and/or Galveston Bay 
water displacement rate, both of which are directly relatable to salinity 
and pollution concentrations. Figure l6 shows a replotted data curve 
from Copeland (l 966) and shrimp productivity by the model for programmed 
variations in freshwater discharge. There is a two year displacement 
(lag) of shrimp productivity change to freshwater input as discussed 
previously. 

Total productivity in the Galveston Bay ecosystem is also a func- 
tion of freshwater inflow. Figure 17 shows the replotted data of 
Armstrong and Hinson (1973) compared to model results for a range of 
freshwater inflows. 

Galveston Bay ecosystem consumer group productivity ratios (of the 
current year to the base-line year) for increased freshwater inflow 
conditions are shown on figure 18. Standing crop biomass levels for 
the six consumer groups are shown in figures 19 through 2k. Table II 
lists the physical and chemical parameters and ecosystem productivity 
as measured by consumer emigrations. 
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Figure 16,- Calibration curve for consumer group 3 due to variation in 
salinity (resulting from changes in freshwater inflow). 
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Figure 17.- Model calibration of total secondary production to freshwater 
discharge to Galveston Bay. 
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Table II: Consumer group and total productivity in Galveston Bay as a result of 
increased fresh water inflow. 


Year 

Freshwater 

Inflow 

cfs* 10" 3 

Waste 
BOD, lb/yr 

*10~ 6 

Salinity, 

PPt. 

Total 

Nitrogen, 

mg/1 

(1969) 1 

IJbTeVoT 

BBSISlli 


1.166+00 

2 

1 . 566 *01 

imfl w 

■ssmsEi 

1.296+00 

3 

1 . 56E +0 l 

2.4 76 * 02 

l. 716*01 

1.296+00 

4 

1 « 566 >01 

2.476+02 

1.716 + 01 

■m-ic.Mi 

5 

1 • 566 *-0 i 

2.4 76 + 02 

1.7 16+01 

1. 296+00 

6 

1.726*01 

2 .476 + 02 

1.656+01 

1.226+00 

7 

1.89E+01 

2.476+02 

1.596*01 

1. 166*00 

8 

2.086*01 

2.476*02 

1.526+01 

1.096+00 

9 

2.296*01 

2-47E+02 

1.466*01 

1.036*00 

' 10 

2.526+01 

2.47E+02 

1.38E+01 

9.636-01 

11 

2.776*01 

2.47E+Q2 

1. 316 + 01 

9.036-01 

12 

3.046*01 

2.476*02 




Year 

Relative Consumer Group Emigration, grams 

1* 

2 

3 

4 

5 

6 

Total 

(1969) i 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

9.33E+05 

5.026+03 

6.876+03 

9.8 76 + 04 


2.42E+03 


1.09E+06 

1.17E+06 

4.61E+03 

4.53E+03 

6.336+03 

5.686+03 

i. 186+05 
l. 216+05 

3.456+04 

4.03E+04 

3.956+04 

3.89E+04 

3.706+04 

Mttlflikd 

1.666+05 
1.74E+05 
1.73E+05 
i. 726+05 
1.626+05 

1.16E+06 

1.16E+06 

1.096+06 

4 .576+03 
4.576+03 
4.706+03 

5-64E+03 

5.356+03 

5.666+03 

l. 216+05 
1 . 2 16+05 
1.126+05 
T.OIE+OT 

9.706+05 

5.056+03 

6.28E+03 

3.2 76 + 04 


1.486+05 

8.52E + 05 

5.286+03 

6.936+03 

Firerarei i 

2.836+04 

2. 486+03 
2.62E+03 
2 .74 E +03 

1.32E+05 

1.166+05 

7. 296+05 

5.506*03 
5. 726+03 
5.92E+03 

6.1 16 + 03 


7.65E + 04 

2.40E+04 

6. 07E+05 

8.61E+03 

6.356+04 

1.956+04 

1. 006+05 

4.916+05 

. 9. 4 26 +03 

5T0 76*04 

1.53E+04 

2.866+03 

8.426+04 

6.956+04 

3.846+05 

wwBfEcmm 


1.146+04 

2. 966+03 









* Number 
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Figure 19.- Zooplankton standing crops during period of increased 
freshwater inflow to Galveston Bay. 
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Model Validation 


Model validation was approached by two independent methods; first, 
a ratio of food consumed to that available in the estuary for each food 
that .is also a .consumer , and second, predicting by the model, effects 
of changes and comparing predicted results to those reported from anal- 
ogous field. studies . 

For the ratio of food (consumer group) biomass consumed to that 
available in the estuary a 10 % biomass conversion factor was used. 
(Pendleton 1973) (e.g. , weight increase of the consumers was 10 % of the 
food consumed). For the base-line year (1969) ratios were 


Food No. /Type 


Food Consumed^ 1 j /Food Available 


U . / Zooplankton 

1.08 * 10° 

5 . /Herbivores 

... 1.13* 10 2 

6. /Omnivores 

1.71 * 10 



7. /Primary Consumer 

1.27 * 10 x 

8. /Middle Consumers 

U,79 * 10 _1 


( 1 ) 


Based on .1 food utilization by consumers (Pendleton 1973) 


Foods U, 5 i and 6 (zooplankton, herbivores and omnivores) repre- 
sent groups of organisms on which we have limited knowledge on growth 
and reproduction rates (zooplankton) and/or incomplete sampling data 
(herbivores and carnivores) since they include small bottom dwelling 
worms, amphipods, ostracods, etc. For foods 7 (primary' carnivores ) 


• . . ; 48 

and 8 (middle carnivores) . the analysis indicates that quantities equal 
to 12. 7 and 47.9 percent, respectively, are consumed in the estuary. 
These values appear to be reasonable based on energy flow- and feeding 
relationships in the estuary, figure 5 (Darnell 1958, 1961; Copeland 
and Fruh, 1970). . 

The second validation procedure, was to predict through the model, 
effects of changes in exogenous variables where results of field studies 
or other analyses could be used for comparison. Two example studies 
on reduced pollution input and reduced freshwater inflow are presented 
in the following section. These are of interest for effective manage- 
ment of estuarine ecosystems. 
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EVALUATION OF EFFECTS OF MANAGEMENT 

The model developed in this study is used to study the effects of 
changes in exogenous variables, specif ically: 

a) decrease in waste discharge to the Galveston Bay from the 

Houston Ship channel (KSC) and 

b) changes in freshwater inflow. 

Waste Discharge 

Armstrong (1973) indicated that due to pollution control measures 
on the Houston Ship Channel (HSC), projected decreases in pollution 
load from this. source were: 

' • • Average 1969-71 • 153*2 x 10 6 lb BOD/yr 

1973 41.6 x 10 6 lb. BOD/yr ; 

Goal (1975) 29.1 x. 10 6 lb BOD/yr 


Effects of this pollution load reduction were evaluated by the 
model based on the following yearly pollution loads: 

Organic Carbon 10^ lb BOD/yr 


Year 

1969 (base-line) 

1970 

1971 

1972 
1973. 

1974 

1975 


Trinity River HSC 
29*9 165.0 

29.9 152.2 

29*9 140.0 

.29*9 80.0 

29*9 41.6 

29*9 33.0 

29.9 


Other 

Total 

52.2 

247.1 

52.2 

235.3 

52.2 

222.1 

52.2 

162.1 

52.2 

123.7 

52.2 

115.1 

52.2 

112.1 


29.0 
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Freshwater inflow was unchanged as were waste discharges from other 
sources. 

Total productivity in the Galveston Bay is projected to increase 
by about bQ% due to the waste discharge decreases, figure 25. This 
effect is due to the overall increased basic foods (detritus and phyto- 
plankton) and decreased pollution toxicity effects on organism's growth. 
As has been discussed previously, total productivity is determined pri- . 
marily by the primary and middle carnivores (consumer groups U & 5s 
respectively); however, the other consumer groups also increase in . 
productivity, figure 26. Omnivores (primarily shrimp) increased rela- 
tively more (6 7$ over base-line year productivity) than the total; due 
to higher overall sensitivity to pollution effects in their own growth 
rates and of their food sources. Yearly productivity values for the six 
consumer groups along with the pollution parameter (average annual 
total nitrogen) are listed in table III. 

Biomass standing crop values in the Galveston Bay during the 7 
year period of the analysis are shown in figures 27 through 32.. As 
would be anticipated from the productivity figures (measured by total 
emigration from the ecosystem) there is a steady increase in biomass . 
standing crop from year to year as pollution input is decreased. 

Freshwater Inflow 

After the pollution load decrease from the Houston Ship Channel, 
which was .assumed to occur from 1969 through 1975 » (see previous exam- 
ple, para. 1), a step decrease of 50$ in freshwater discharge from the 
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1969-71 Average 
(Armstrong and Hinson, 1973) 




Figure 25.- Projected HSC pollution load and effect on Galveston Bay 
productivity. 
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Figure 26.- Galveston Bay productivity ratios for six consumer groups 
and total "biomass due to decreased waste discharge to the 
Houston Ship Channel . . 


Table III: Consumer group and total productivity in Galveston Bay as a result of , 
decreased waste discharge to the Houston Ship Channel. 


• 

Relative Consumer Group Emigration, grams 


1* 

2 

3 

4 

5 

6 

Total 

(1969) i 

9.336+C5 

S.f>?6*0 3 

6. 8 76 + 03 

9.P7E+C-4 

3. 086*04 

2.4?6*C3 - 

1.446 + C6 

2 

9.475+05 

5.1 If 40 3 

* .996*07 

l .0 1 6 +C 5 

* . 1 2 6+04. 

2.44 6*0 3 

1 .465*05 

3 

9.7of*C8. 

5.256+03 

7. 716*03 

1. 0 7 C + 0S 

3. 22E+04 

2.53E+C3 

1. 51E+C5 

4 

• 1 .06 C +Cft 

5. R2F + 03 

3 .?2 C *03 

1.155+05 

3 .506*04 

2 .756*03 

1 .676*05 

5 

l .18 C +C6 

*.51503 

9.9 76*0 3 

t . 78F+C5 

4. C8 C + 04 

3. 246+07 

1 .896+05 

6 

1 .2'5frf*6 

5.83E+03 

I .1 1 c *04 

l .34F+05 

4 .44F+04 

3 .5 R 6*03 

2.005+05 

7 

1 • ? 8 F +0 6 

6.95 5+C3 

1 . 146*04 

l. 376+05 

4. 516*04 

3. 7QP + C3 

2.045+05 

8 

1.28E+06. 

6.9*5*03 

1. 146*04 

t .3 7 c +C5 

4.536*04 

3 .746+03 

2 .0 45 + r 5 

9 

1.786+06 

6 . 9 5 c +o 3 

L . 14 p+04 

1.376+05 

4.616+04 

3. 736+03 

2. 046 + C 5 

10 

l. ?R c +06 

6.956*0-* 

1.146*04 

1 .-37£*OS 

4.516+04 

3 . 73 c +0 3 

2,.C4F-+j6 

11 

1 ,2« c +0 6 

5 .95 50 3 

1 . 146*04 

1.375+05 

4 . 6 1 c +04 

3.736+03 

2.046+^6 

12 

t. ?fl c +CS 

6. 9 5 C + ■- * 

i . 1 4 C *04 

l . 3 7F f 0 5 

4.51E+C4 

3 . 73 c *0 3 

; — 

2.046 K 5 


Year 

Freshwater 

Inflow 

1000* s cfs 

Waste 
BOD, Ib/yr 

♦ 10 -6 

Salinity, 

ppt. 

Total 

Nitrogen, 

‘ mg/1 

(1969) 1 

1 .096+7 i 

2,475+07 

I.59E+ 01 

i. 1*F*0P 

2 

1.895+01 

2.35f+^2 

1.596+01 

T.lTF+00 

3 

1 . oqf +/>1 

' 7.77F +C’ 

1 . 896+01 

. 1. 1CF+00 

4 

1 .895+^1 

1 .*28* )? 

1 . 595+01 

Q.36F-01 

5 

1.895+01 

1.P4F+0? 

1.SQ5+01 

7.R3F-01 

6 

1 . P 9 f +0 1 

1 .1 5 = 07 

1.595+0! 

7. 78P-0] 

7 

1.856 + 01 

t. 126 + 07 

1. 59C+01 

*.906-0] 

8 

1 .P9F+01 

1 .1 7 5 + 9 7 

1.596+01 

6.996-01 

9 

t.por +71 

1.175+07 

1. 595+ 31 

6 • 9OF-0.1 

10 

1 .P95 + M 

1 • ! 2 c + 9 2 

1 .89p*C 1 

6.99F-01 

11 

1.896 + 01 

1 .12 P *C2 

1 . 80 5 + 0 1 

6.996-01 

12 

1.898+4-1 

1 .126*07 

1.895+91 

6.996-01 



♦Number 
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Figure 31.- Middle carnivore standing crops during period of reduced 
waste discharge to HSG. 
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Figure 32,- Top carnivore standing crops, during period of reduced 
waste discharge to HSC. 
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Trinity River (e.g., due to filling a reservoir) was assumed for the 
calendar year 1976, followed by a restoration to 75$ for the calendar 
year 1977 and three following years. Freshwater inputs to the Galveston 
Bay were then: 

Freshwater Inflow, 1000 CFS 


Year 

Trinity River 

HSC 

Other 

Total 

1969-1975 

5.58 

10.53 

2.79 

18.90 

1976 

2.79 

• 10.53 

2.79 

16.11 

1977 

h.lQ 

10.53 

2.79 

17.50 

1978-1980 

h.lQ 

10.53 

2.79 

17.50 


Total, productivity in the Galveston Bay ecosystem increased about 
30 % above the base-line year due to the combined effects of reduced 
waste discharge to the Houston Ship Channel, followed by reduced 
Trinity River freshwater discharge, figure 33 and table IV. Effects 
due to the first change (reduced waste discharge) are discussed in the' 
previous section. Reduced freshwater flow leads to increased total 
productivity; however, one of the commercially important , consumer groups 
(omnivores, which includes the shrimp species) as reduced about 20 %, 
due to. their sensitivity to system changes (salinity and pollution con- 
centrations). The middle carnivores (consumer group 5 which includes 
the anchovy) are favorably affected (increased about 1&%) due to 
decreased freshwater and the consequent . increase in salinity. The 
other consumer group productivities are moderately increased (about 10$) 
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Decreased HSC waste 



Figure 33 «- Effect on Galveston Bay productivity due to Decreased HSC 
Pollution load (1969-1975) followed' ty decreased Trinity 
River freshwater discharge (19.7-6- 19 80) • 


Table IV: Consumer group and total productivity in Galveston Bay as a result of 
decreased waste discharge (in the HSC) followed by decreased fresh- 
water flow from the Trinity River. 


Year 

Freshwater 

Inflow 

cfs* 10" 3 

Waste 
BOD, lb/yr 

♦ io - 6 

Salinity, 

ppt. 

Total 

Nitrogen, 

mg/1 

(1969) 1 

— i>d9Et01 

__2.4 7E402 

1- 59F401 

! . 1 6 F 4 R 0 

2 

.„U£3£».flLL-j 



1.1 3F400 

3 

..1-j89_£> 01 j 

Z.ZZLtD^ 

1 . 59 F 401 


4 

_Ud*j=*o±_ 

1 . t>2f *-02 

i . 59F401 

9.36F-01 

5 


A .2 4 £*-02 

1 . 59F 401 


6 

UQ2L* 31_ 

1 . 15E *02 

1. 59F*01 

7 . 2 8E-0 1 

7 

U*fl9E*Ql 

1.I2E40? 

1. 59E*01 

6 . 99E -01 

8 

- 

-L.J.2E 4 3 2 

L. TOE * 01 


9 

1*_75£ 4QJ_ 

~-L 


7. 28F-01 

10 

- 1*..2££±j0lL 

-JL .12E_ta2_ 


7 . 2 8F-GI 

11 

-i . -tUL. 

.j..Z£±nz_- 

4Q1 


12 

01- 



—L» 28 E^LL_ 


Year 


Relative Consumer Group Emigration, grams 

1* 

2 

3 

4 

5 

6 

Total 

(1969) 1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

9.53E4 J5 
9.47E405 

5.028*03 

-S*X1E*J3 

£^£.*03 

9.87E 404 

3 • 08E +04 

2. 42 F 433 

1 .4.46*0 5 

6.99F403 

1 .01 E 43 5 

3 . 1 2E * 04 

2.44E403 

1 .46E405 
1 .51E405 
1 .678405 
1 . 89E 40 5 

9. TOE 4 05 

5 «25E 403 

7 . 3 1 E 403 

1 • 03c 4 05 

3.22E404 

2 . 53 £ 403 

1 .06 E 406 

5.82L403 

8 . 22E 4 33 

l .15E*u5 

3.50E 404 

? . 75F 403 

1. 186406 

_6 ^5 IE *03 

.9jj_97E. 4 03 

1.1 IE 4 34 

1.28E405 

4 • 08E 434 

3 . 24F 403 

1.25 £436 

j6«J3E»Q3, 

l .34E435 

4.446404 


2 . 006405 

1 • 28E 4 06 

6 . 95JE 4.0.3 

1..J.4E * 04 

T.37E405 

4 • 5 1 E 404 

3 .7 OF 403 

2 .04E 405 

1 .49F406 

6. 73F 403 

1.1 IE *04 


5.02E4 34 

3.68E4J3 
3.63E403 
3^,68 E *03 

2.326405 
2 .266405 

1 . 4 6 F 4 06 

6. 78F433 
6.86E403 
.^9_4E_4 03 
6 . 84E 4 03 

_1 .Q8F *04 

L. 51 E 4 05 

5.38E404 

1.408*06 

1 .49E405 

5 .00E 404 

2.206405 

2.206405 

2.2Q£»05 

1 - 4 IE 4 0 6 

1 . 07 F 4Q4 

1 . 49E * 05 

4 • 95 E 434 

-LmJOA+SA 

I » 1 OE ♦ 04 

l .49 E 40 5 

4 .97 E 4 04 

3.67F403 

1 







♦Number 


o\ 

ro 
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by the decreased flow. Productivity ratios (ratios of current year to 
baseline year) for the six consumer groups and average annual salinity 
and total nitrogen concentrations are shown in figure 34. It is inter- 
esting to note the carry-over effects in. the ecosystem predicted by the 
model. For example, there is a strong perturbation on consumer group 
3 (omnivores) due to effects from prior years. The specific reason for 
this lag is unknown and may be due to reproduction or food effects , 
but has been well documented in previous studies (Copeland 1966).. The. 
significant idea however is that short term (year to year) comparisons 
may be misleading where management of an ecosystem is concerned. 

From an estuarine fisheries viewpoint, biomass levels within the 
Bay are important since they represent concentrations available for . 
harvest' (catching) at a given time. Biomass levels for each of the six 
consumer groups 6ver the period for the assumed conditions are shown in 
figures 35 through 40. Zooplankton density increased in 1976 and 1977 
due to decreased freshwater inflow. Subsequently, short term effects 
of increased foods and increased pollution appear to offset each other, 
with a subsequent stabilizing at a lower level. As indicated from the • 
productivity analysis the omnivores (consumer group 3) are the ; most 
sensitive to system changes. There is a sharp decrease in biomass 
levels of these consumers followed by relatively rapid .. stabilization 
at a value about 10# below the point before freshwater inflow was 
decreased. Consumer group biomass levels tend to shift about due to' 
delayed responses in growth of foods within the estuary,, then stabilize 
at the new level after about 2 years. One of the factors that must be 
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taken into consideration, of* course, is that year to year effects are 
not completely understood and thus possible effects have been minimized 
in the’ model due to 'lack of quantitative data. 
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DISCUSSION AND RECOMMENDATIONS 

A simulation type mathematical model that quantitatively relates 
the chemical, physical and biological characteristics of an ecosystem 
has been developed for the Galveston Bay, Texas. The model incorporates 
results of a number of investigations into detailed processes in the 
ecosystem including those of seasonal immigration and emigration, food 
and consumer densities , feeding habits and responses of organisms to 
exogenous changes such as waste discharge and freshwater inflow. 

Two examples of possible management actions are analyzed using 
the model. Comparison of model predictions and results from analagous 
field studies over a 20 year period and reported in the literature demon 
strate the model's usefulness. The first example evaluates the effects 
of reduced pollution inflow from the Houston Ship Channel. Increased 
overall productivity with relatively higher increases for shrimp is in 
agreement with Copeland (1966), Armstrong and Hinson (1973) and 
Wohlschlag (1972). The second example investigated the effects of 
reduced freshwater inflow. Predicted results of increased zooplankton 
and finfish and decreased shrimp are in agreement with Copeland (1966), 
Cooper (1967), and Armstrong and Hinson (1973). Additional advantages 
of the model, of course, are quantification of the results and more 
complete analysis of time-dependent (system lag) effects. 

One of the purposes of the effort was to develop a model that may 
be used to provide predictive information to those responsible for 
estuarine management. To the extent that a quantitative relationship 
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has been developed and verified (at least with historical data) this 
basic objective has been accomplished. The degree that those quantita- 
tive relationships are used will be determined by the confidence and 
subsequent use for inputs to management action. Obviously, determining 
the degree to which this has been successfully accomplished is beyond 
the scope of this particular effort. 

One factor apparent from this modeling effort of ecosystems is the 
general lack of in-depth studies of the multi-parameter effects of 
seasonality, food , temperature , pollution , environment , etc. on the 
survival and growth of estuarine species. Laboratory studies have been 
made in some cases of the better known and/or commercially important 
species; however, due to various limitations these have seldom been 
extended to the ecosystem level and when they have, with limited con- 
sideration of more than one variable; the others held constant , ignored, 
or assumed of secondary significance. 

In this study, only those variables subject to control and conse- 
quently important in management decisions are classed as primary vari- 
ables (e.g., freshwater inflow, waste-discharge and their distribution 
in the ecosystem). It is probably important to know the effects of 
such variables as temperature (to assess power plant locations, for ' 
example) but its significance as a management tool may be limited, 
particularly if its control is beyond the scope of management. (On the 
other hand, knowledge of these effects may be necessary to develop a 
model with adequate predictive accuracy for power plant siting studies, 
etc. ) 



In the present study emphasis has been placed on ecosystem response 
as opposed to detailed responses within the ecosystem (eig. , of specific 
species). Thus, foods and consumers were grouped on the bases of func- 
tion rather than studies of individual species of plants and/or orga- 
nisms. This has many advantages from the modeling viewpoint s . but ob- 
viously leads to problems for the biologist, since his groupings may 
not be the same and/or there are interactions within or outside of the, . 
group that are not compatible with laboratory investigations. On the 
other hand, it is obviously necessary to know more about the ecosystem 
than relationships such as ... "when freshwater input increases , shrimp 
catches increase two years later". The missing information -is a quan- 
titative description of the links in the energy flow over the two year 
period. For example, what ecosystem processes .are affected during the 
delay? Are these due to hatching and spawning, immigration, food avail- 
ability, predator-prey relationships, etc? In addition, each of these , 
quantitative relations must be compatible with interlocking processes in 
the ecosystem. Thus, in theory, an analytical description of an eco- 
system is an almost infinite number of simultaneous equations.' 

With some background and experience in assessment of ecosystem 
response it is possible to quantify some of the more significant rela- 
tionships. This is. the basis of this modeling effort. From additional 
studies of the model the critical parameters may be determined and these 
should be the basis for further studies to improve our understanding and 
predictive capability of the model. - 
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As a result of this modeling effort which uses analyses and data 
from previous studies, additional data are required in the following 
areas : \ • 

1. Affect of pollutants on survival and growth of dominant estua- 
rine species ; 

2. Absolute or relative effects’ of pollutants on different' food’ 
and consumer species in an estuarine ecosystem; 

3. The degree and extent of variations in immigration, and emigra- 

from year to year; and ! 

h. Quantitative definition of observed biological, time delays 
1 such as the approximate two year lag in shrimp productivity 
•after changes in freshwater inflow. 

, The list is only limited by the accuracy requirements of the modeling 
processes; that 'is, the objectives of the program. . 

It is recommended that the simulation model which includes analyt- 
ical and empirical descriptions of Galveston Bay ecosystem processes 
be used where possible to predict effects due to exogenous changes in 
-freshwater inflow, waste discharges, and other system parameters. When 
results of later studies become available these should be incorporated. 
Further, the modeling and ecosystem processes concepts developed in 
this study should be applied to other ecosystems, whether they are 
subject to natural or man-influenced exogenous effects. 
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CONCLUSIONS 

As. a result of this investigation and analysis of the results the 
following conclusions are made: 

1. A simulation type mathematical model has been developed, for the 
Galveston Bay, Texas, ecosystem that quantitatively relates . 

t . • . . 

pollution and freshwater inflow to secondary productivity. ' 

2. Ecosystem responses for two potential management actions, 
reduced pollution input and .changes in freshwater inflow, were 
evaluated and results may be studied by management to determine 
the desireability of such actions. . 

3. Outputs of the model are quantitatively applicable, to the 
Galveston Bay. In the same temperate zone (approximately same 
latitude and weather conditions ) the model is readily adaptable 
to other ecosystems. Initial numerical results should be 
reviewed carefully for local effects. 

4. For other temperate zones and/or different environmental (rain- 
fall, tides, etc.) conditions a study must be made to determine 
similarity among consumers, food types and availability, and 
seasonal growth characteristics between the Galveston Bay 

and the ecosystem of interest. Calibration and model verifi- 
cation will probably require a comprehensive sampling effort 

in addition to historical systems responses (sampling or 
commercial records in conjunction with rainfall, for 

example ) . 
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Development of a reliable predictive model is an iterative 
process that improves with iterations; the systems analysis 
model provides an optimized framework for analysis using all 
available information. 
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APPENDICES 



APPENDIX A 

Functional characteristics of the Galveston Bay model are described 
in the model section. Consumer group biomass values are developed in 
two basic steps. First, productivity (biomass) curves for the baseline 
year (1969) use well-established ecosystem characteristics (£*£.* , 
immigration, growth and emigration) to determine, by iterative techni- 
ques, model constants. The basic objective is reasonable agreement 
between the model biomass curve and sampled data, figures Ik and 15. 
Second, for years after the baseline, biomass growth rate is the pro- 
duct of the prior year's net' growth rate ’and growth rate ratios (of * 
the current to prior year) due to changes in 1 

1) food and consumer densities; 

2) environmental factors (including pollution); and 

3) distribution. 

Equations explicitly describing the above functions are discussed in . 
the following sections. ‘ 

Baseline year 

Baseline year biomass curves explicitly describe sample values. arid 
biological factors in model language. The pertinent ecosystem factors 
(e.£. , immigration, growth and emigration) are shown schematically in 
Figure Ik. Reference will be made to curves "a" through "e". Model 
baseline year constants for consumer group b are listed in the model 
section. Model equations are expressed in terms of the calendar year 
periods, I. 
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Immigration (curve "a") is taken as the positive portion of a sine 
function ( tt radians) and is described by the equation 

C0NSIM(I,L) = SIN(R1(L) * 3 .1U16/R2 (L) )* CONSTIM (L) (l) 

where 

Rl(L) = number of periods since start of immigration! 

R2(L) = number of periods of immigration; and 
CONSTIM (L) = maximum immigration level, grams 
(where Rl(L) > R2(L), CONSIM (l,L) =0.) 

Larval and post-larval populations have very high growth rates, 
decreasing with increasing organism size (Patten, 1971 i Paloheimo and 
Dickie, 1965). A decreasing, value exponential function is used to • ' 
describe the growth rate of immigrated consumer groups in the baseline 
year (adult biomass in the baseline year has a growth rate of 0). The 
model equation is 

•EXPGR (M,L) = ( [ exp ( . 05 * (27 - M) ) ] - X.) * BLYGRL(L) (2) 
where 

BLYGR * the growth factor for consumer group L in the baseline 
year. 

#GEPYR(I,L,1) in the model when correlated with the calendar 
year periods. 

Period growth rates are determined by this relationship for each 
consumer group to obtain the baseline year biomass curves . Growth rate 
multiplier is 1.0 plus growth rate. 

■ . . 

After growth in the estuary the consumer organisms emigrate from the 
ecosystem - in this case primarily to the Gulf of Mexico . In the 
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model emigration is the sura of two components - early and late emigration. 

a. Early gration EAREM(I,L) - is that phase of emigration which occurs 

at about the time of the biomass peak of that consumer group and is 
determined by: . 

*EAREM(T,L) = SIN(R3<L) 3.lUl6/R4(L) * FACTEM(L) * YCONSIL(L) 

' • . ' ... . - (3) 

where . 

R3(L) = No. periods since start of early emigration; 

RU(L) = No. periods over which early emigration occurs; 

FACTEM(L) = Constant multiplier for consumer group; and 
YCONSIL(L) = Biomass of young of consumer group L, grams, 
for the period of calculation. 

(IF R3(L) > nh(L). EAREM(I,L) =0.) 

b. Late emigration (F0REM(I,L) is determined by 

#F0REM(ljL) = YCONSIL(L) * (R5(L)/ (R6(L) + 1 . ) ) (4) 

R5(L) •* Periods since start of late emigration - 
R6(L) - Periods from start of late emigration to end of 
seasonal year. 

(IF R5(M,L) > RefMjL) F0REM(I,L) = 0) 

•In the model the total emigration is C0NSEM(l,L) 

Numerical values for the constants in the above equations were determined 
iteratively to provide reasonable fit between model biomass curves and 
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sampled data. An important step in the baseline year portion of the 
model is to sequentially correlate the individual consumer group sea- 
sonal periods (M) with the calendar year periods (I). Note that in the 
baseline year portion of the model the consumer group biomass curves 
are developed independently. 

In the model, consumer growth rate multipliers, and ratios are • 
designated by the acronym GEPYR (I,L,K) for young » this years 

immigration) or GEPAR(l,L,K) for adult biomass where 
. ! _ ' 

I = calendar year model periods, I = 1 ,..., 2 6 ; 

L = consumer group, L - 1,... ,6; and 
K = sequential number for growth multipliers 
K =1, for baseline or prior year" 

K = 2, ratio for food and consumer densities; 

K = 3, ratio for environment (including pollution); 

K » 4, ratio for distribution'; and 

K = 5, net growth rate current year (product of 1-4 above). 

GEPYR (I,L,l) values are determined for the baseline year using values 
from the exponential curve (In the baseline year all other GEPYR and 
GEPAR multipliers are equal 1.0). Thus GEPYR(l,L,5) equals GEPYR(l,L,l) 
and GEPAR(I,L,5) equals 1.0 for the baseline year. In subsequent years 
GEPYR(I,L,l) for the current year is set equal GEPYR(I,L,5) for the 
immediately prior year; GEPYR(I,L,K) , K = 2,..., 4; values are 
calculated; and GEPYR(l,L,5)for the current year is the product of 
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GEPYR(I,L,K), K = Adult growth rate multipliers are deter- 

mined in an analogous manner . 

Predictive Model 

After the baseline year, biomass levels in the ecosystem may change 
due to changes in net growth rates, which result from new values of food 
and consumer densities and/or exogenous' variables . Revised net growth 
rates (the product of prior year net growth rate and growth rate change 
ratios) are used to determine current year biomass levels by iterative 
calculations. . 

Exogenous variables (inputs) in the model are freshwater inflow and 
pollution load. The former directly affects only salinity and both 
variables affect the pollution indicator, total nitrogen Oe. , by dilu- 
tion or pollution load) . Average annual values are used in the model as 
discussed in Section II. These physical and chemical changes are used 
to determine revised growth rates due to 

1) Food and consumer densities; 

2) Environmental factors (including pollution); and 

3) Distribution, 

which will be discussed in that order. 

Food and consumer densities . 

Food and consumer density effects on growth were studied by Brocksen, 
Davis and Warren (1970) in an environment where there was both a dis- 
tinct consumer and a distinct food. However, in the Galveston Bay 
ecosystem there are a variety of consumers and a number of foods as 
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discussed in the ecosystem section. From the feeding preferences of the 
Galveston Bay consumers energy source matrices were developed from 
iigures 6 through 11. This figure (^l) shows a typical period of a 
matrix which is 2 6 ( e . g . , I) segments deep. There are separate matrices 
for the young and adult consumers. Note that the columns add to 1.0 since 
they represent the total energy intake of each consumer group. The 
purpose is to provide a model-period relationship between consumer food 
requirements (or desires) and food density. Further, by comparing the., 
current year's consumer group biomass to the prior years, we may assess 
consumer density effects. 

Explicit expressions used to, determine food availability, food and 
consumer densities and resulting growth rate ratios follow (Note that 
ic.ods i through 3 are not affected by consumer group growth rates, but 
foocs 1* through 8 are consumer groups 1 through 5, respectively). 

Armstrong & Hinson (1973) estimated that 5&% of the detritus 
organic carbon comes from marsh and other rooted grasses,. In this 
study it is assumed that the above grasses are influenced by the same 
factors that affect the growth rates of phytoplankton ana detritus • 
if from the prior year’s growth (Heald, 1973). Freshwater flow, rate 
from the rivers will affect the flushing rate from the marshes and 
rivers and is accounted for by a .5 power factor on. the water flow 

f * 

ratio; which is applied to the remaining h2% of the organic carbon. 

There are no prior studies in this latter area; however this is not a 
sensitive factor in model performance. 

Detritus for years after the baseline year is determined by; 
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Figure 4l.- lypical period of Galveston Bay consumer energy source 
matrix. 
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FOOD(N,I,l) = FOOD(K,I,l) * {.k2 * (DWF ** .5) + *58 * PFHG) (5) 
where 

F00D(N,I,1) = detritus {food l) for the Ith period of the Nth 
year; BOD, lb/day , 

F00D(K,I,l) =' same as above except K = N - 1; •• 

DWF = ratio of total freshwater inflow, current (N) to prior 
(N-l) years; . ‘ 

PPGR = ratio of phytoplankton growth rates of prior (N-l year) 
to two years ago (N - 2 year). 

Phytoplankton concentrations are a function of short term (seasonal) 
environmental parameters! Growth rates are based on data of Copeland • 
and Fruh (1970) from which the following relation was obtained by a 
multiple-regression analysis: 

PPG(N) = .8211- .0207 * AANIT(N) - .0129 * AANIT(N) 2 + (6) 

.0l*l9 * AASAL(N) 

where 

PPG(N) - Phytoplankton growth multiplier in Nth year 
AANIT(N) - Average annual total nitrogen in Nth year, mg/L 
AASAL(N) - Average annual salinity in Nth year,ppt 

Vascular plant material - due to the large quantities of submerged 
grasses in the Galveston Bay this material is always considered to be 
in excess. 
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Foods 4 through 8 are consumers in the Galveston Bay and current 
biomass levels are used in the iterative calculations. Food 4 quantity 
equals biomass for consumer group 1, food 5 for consumer group 2, etc. 

Only "young” consumer levels are used in determining the equivalent 

i 

food densities. . 

Food density effects are determined' from the investigation of 
Brocksen, Davis and Warren (1970 ) (see Figure 13) and are defined in 
the model as follows : 

_ 1 Hr .01 » (l.2».75»A LOG 10(FOOP(N,IH,J)»FQ(J))) / ?) 

FFWJ ~ ± + .ol#(l.2*.75*A LOG 10 (FOOD(N-l) ,IH,J)*FQ( J) ) ) 

where 

FF(J) - growth multipliers due to food J; 

F00D(N,IH,J) - Food quantity - Nth year, IH (= I - l) period, 
Jth Food, BOD (lb/day) or grams; 

F00D(N-1,IH,J) - Same as above for prior year and 

FQ(J) - equalizing factor for various food densities and types. 

Specifically, constants FQ(J) in the expression are related to maximum 
food density of U00 mg/m^; thus FQ(j) values are 1+00/ (Max level of food 
J in baseline year) (Brocken, Davis and Warren, 1970). 
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The term FF(J) is then multiplied by the appropriate energy source 
matrix (ENERS) term and summed over J = for each consumer group. 

This food-effects factor is designated R(L) in the model. 

Quantitative expression for consumer density is (Brocksen, Davis 
and Warren, 1970 ) : 


GRCONSD ( L ) = 


1 + .0215 (1. - (C0NSIL(l)/C0NSTL(7)» 4.) 
1 + .0275 (1. - .25) 


( 8 ) 


where 

COKSIL(L) is the consumer biomass level at the beginning of the 
current .period, grams - 

C0NSTL(7) is the consumer biomass level for the same period in 
the base-line year (at biomass levels four times those of the 
baseline year, growth rates are taken as zero; Brocksen, Davis 
and Warren, 1970), grams. 

Calculation procedure is to use food and consumer densities from 
the previous period along with the energy source matrix for the current 
period (note that food densities at the end of this previous period are 
the same as the values at the start of the current period). 

To allow for the omnivarious feeding of estuarine species (Darnell, 
1958), if the food - consumer density parameter is less than 1, the 
calculated value is raised to the .8 power (i.e., hungry organisms make 



more effort and/or are less food selective); no. quantitative- work is 
reported; however, the model is not sensitive to this value 

The food-consumer density parameter (product of R(L) and GRCONSD(L)) 
is raised to the 1 if— power, which is the number of days in a model 
period. This term is designated GEPYR(l,L,2) (or GEPAR for adult bio- 
mass) and is the growth; ratio compared to the prior year. , 

Growth change for the six consumers groups due to environmental 
parameters changes is 

GEPYR ( I , L , 3 ) = FG(N)/FG(N - l) (9) 

where 

GEPYR(I,L,3) - growth ratio of current to prior year due to . 
environmental parameters change (young organisms); 

FG(N) - fish growth multiplier in current year (FG(N) = 1.1 - .1 *{l0) 
AANIT(N)/AANIT(l) ; where AANIT is average annual total nitrogen) ; and 
FG(N-l) - fish growth in prior year. 

Growth rate effects due to changes in distribution of organisms in 
the estuary as a result of changes in environmental parameters and their 
relation to geophysical (and consequent water and waste distribution) 
factors is important in the life and growth of estuarine organisms 
(Copeland 1966, Armstrong and Hinson 1973)* Temperatures and other 
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factors are apparently, of secondary importance relative to salinity and 
pollution effects in this ecosystem in successive years. In the model, 
a quantitative relation is based on the concept that when freshwater 
input decreases shrimp and crabs (consumer group 3) have less area in 
which to spawn and grow with minimum predator effects i while zooplank- 
ton and primary and secondary canivores -have greater feeding area and 
a predator advantage due to greater penetration into the marshes. 

Empirical equations are developed from data and analyses of Copeland 
(1966) and Armstrong and Hinson (1973). 

' .0 . 

For omnivores (consumer group 3, shrimp and crabs) 

GEPAR (1,3,4) = (1. +( .54*( (OMNG(N)/OMNG(N-l))-l.),)‘) ## (l./26.) (U) 

where 

GEPAR(I,3,4) a growth rate ratio in period X for consumer 
group 3 due to growth changes of this consumer group due to 
distribution factors (sequential number 4); 

OMNG(N) = omnivore (consumer group 3) growth rate for the 
current year and is determined by OMNG(N) = 0.09467* 
(25.-AASAL(N-2) )**2. 05 (see Figure l6); and OMNG(N-l) is 
same as above for prior year. 

Galveston Bay total productivity (biomass) is dominated by the 
consumer groups 4 (primary carnivores) and 5 (middle carnivores). ; 
Zooplankton aiid herbivore growth rates would be increased as a result 
of increased salinities (Cooper, 1967). Growth rate ratio equations are 
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GEPYB ( 1 , 1 , 1* ) = .50 + .50*AASAL(N)/AASAL(l) (.12) 

where 

GEPYR(l»l,^) = growth rate multiplier in the I— period for the 
first consumer group due to distribution 

effects (h‘) and • 

AASAL(N) - average annual salinity, current year 
(Note that average annual salinity for baseline year is 
AASAL(l) . ) 

GEPAR(l,l,U) , adults of zooplankton; GEPYR(i,l+,U) , GEPAR(I ,U ,U) , young 
and adults of primary carnivores; and GEPYR(l , 5, J 0 , GEPAR(I,5,^) young 
and adults of middle carnivores are determined in the same manner . 
Application of the above empirical equations are shown in the calibration 
curve. Figure 17. 

A net growth rate is determined from the product of the prior year . 
growth rate and change ratios due to 

(a) food and consumer densities; 

(b) environmental parameters and 

(c) distribution 1 

This product is determined and applied iteratively period by period in 
the model to obtain new biomass levels for the six consumer groups. 

Adult biomass levels for each consumer group are based on the same 
factors as for young, except using energy source matrices and food 
curves for this particular class. In the baseline year, it is assumed 
that; due to mortality, exchange through the passes connecting the 
Galveston Bay to the Gulf of Mexico and other factors, the net growth 
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rate is a constant value, equal to 0. (i.e,, multiplying factor = 1.0)- 

In subsequent years, emigration of adult biomass out of the estuary is 
assumed to increase in direct proportion to the ratio of net growth 
rates of the current to previous years (jLe. » if too competive for 
food, adults leave). . Biomass grovrth greater than the change in the eco- 
system emigrates and is included in the total emigration. 

In the model, adult biomass change in a period is determined by: 


PERGRA(I,L) = ABGROW * (GEPAR(l ,L,5 ) - 1.0) (13) 

where 

PERGRA(l,L) - adult biomass change in period X for consumer 
group L grams; . 

ABGROW - adult biomass at beginning of period I for consumer 
group L grams; and 

GEFAR(l,L,5) ” net growth rate multiplier 
Adult biomass at the end of the period is 

ACONSIL(L) = ABGROW GEPAR(l,L,5)/FEPAR(I,L,5) (1*0 

where 

ACONSIL(L) = adult biomass grams at the end of period I; and 
FEPAR(l,L,5) - net growth rate multiplier for prior year. 

Adult biomass migrating during the period is 
PERAFM(I,L) = PERGRA(I,L) - {ACONSIL(L) - ABGROW) 
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This latter value is added to the "young” biomass emigration to obtain 
the total emigration for the period. 
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APPENDIX B 

A computer program print-out, including input data and coefficients, 
is listed for the case of increased freshwater inflow to the Galveston 
Bay. Results are listed in Table II and plotted in figures 18 
through 2k . 
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JMfcNsjUN ExP(pM 2S*6) OOM? *6W (b) *KN<?fe*6) 

J 1 A E M s i ( » \J t K ‘ l ( ^ 1 ) * t '-12 1 O • r.' 1 • 1 ; M i :■» •' ei ) * P ( 2 S • b ) » F OU U I 1 /?(» • H ) 

YMSOMM 

OMENS JON t'-Atf*T <26«6) ♦ CONSUL ( 26 •'!>>-« CON ^S.( 26 • f> ) *C0M^6 

1lMS»*>) (2‘)M6*s) * F d'YO ( 20 • 2S ♦ •* . ) 

Ji Me NS TON rtfilf y^U (20***) .•*“>TtDU (20 *4) MWT ( 20 * 4 ) MPrfT ( ?0 * 4 1 MWP { ?G • 
14) .*VAPAM( 7,27) 

>1 MENS ION vs M6M,r> > *StC0N».(26*h) *w6 <S) *FJ <M) 
j Mr: NS I UN *>> ) I SU ( 20 ) , *Tp]s H2U) *#'»>.il ( ?7 ) • WTQH «?7 ) 
ji -ttNsroN in (.2) * yc:*mp<S20> mp<s?o> , ycjnpj (-> 20.6) 

) I Me. NS I UN O'* NO (2-0) MP(2h) 

) Mr. Nil ON IM6)*I2l6>«I.Mt>>*l4l6>.isl6>«lM6),-<Me»),W?f6),-M(6>,K<f 
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2> • CO OU(^.6hC0NM“(f , r.F) * LJNPlM ( 2 b » S ) *FpP£M ( 2b*6) * TONS M ( b ) • TC 
3.NNSLM ( h ) ,C.) h ibT Im( 6) ,FACTF M (s) , SAMP.(2 cm h ) *.(VFPA* (2b«o*6) « Gt* 3 YP ( , b , 

46) ML Yrw <6) . PC ALP AM 1 dh . t> ) •COM* (M « f AV ( H ) ,-F F ( tf ) • AP <6 ) , AM ( 6) 

->• FKP Yw (26*6* 6) ,K£P Aw 126,6 «S) MP<>(2U) « Ku ( 20 ) • T F *0 1 2 0 ) ,T*D<?U> . AASAl 
5 120 ) 

6 , A AN I T 120) MC0NpFM( 6> • wt SJLT < ?0 • 6*'?01 *UAT A ( 20 • I 03 ♦ PE WG* A ( 26 • *U . 
7- 3 FN, j -,y (26*6) ,T4FC<>NSM) , AVDF JlDM) * TPF CONS ( 2r> , H ) « PFCONP < 2b • « • 6 ) . 
•l^CAUl .ACJOSOMlb) 

j l m r N S I u n C0NS7LM0MS.S) «-»rVAPM(P6,6> ♦ Y m { ? ) 

CALL wbtufiO 

INI 1 > = 1i)H 12 i2 PW.JH - • : • 

M (2) =10**206 PA TO 

JPsii 

XMsJfjUYFflP 
=»FAF = rt.oS. 

J FYF = 0,'.is 

neau cs.i«n (F.)i) 0 (i,ui),i = u?h) 

- 10 FQw-Mat (MFlO.O) 

JO 262 MIMS \ 

*£ao ISMS 1>.f D on ( 1 . 1 ,2) 

263 FOP “1 at (+10.0) 

F JOJ 1 1 • r m m i . 

^s? con r i Mot 

: JFAO In PA-lPLfJ vAL’JFp 
: SAMPLE . ) VALUE p In 3ALV SAY s(I*L) 

JO s.H. Ml MS 

PEAn 1SM0)S{ I , 1) ,b< 1M) f S( M 3) «b( r*4) •su ,S) Ml Mb) 
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FO**AT (6F10.U ) 

5UiM=><I*6>/ia.. 

IFU.GT.DOO TO bO 
«>RJTt <6.61 ) 

61 FORMAT ( 1 ^xS'loAMfVbX-j'-oAM^ 3 b* bMbfl ^4S X 5M5 AMHbbX =>Hb AM^h > 

60 COMT iN'Jt 

#^lTt-{6*6<^) d ( I ♦ 1 ) « S f T . £ ) . b < I . - 3 ) * b ( I * 4 ) » b ( I » S-) * > t I . 6 ) 

Kl Fl>MflT U X *6f 10.0 > 

brt CO'JT I m 'JF 

c fcNEQM** J.U M otuithj aFTF-' ! * • J "OOD.L CONb jM t ;R . 

00 67 *=1 *^6 .'. 

00 o7 J- 1 ,6 
)0 ->? L=1.6 

J.L) =0. 

S7 COMTI'iJL 

30 ^4t) ] t Pb 

^M=FlOAT (m) 

1 , 1) = . 1 
tMt -<=> < •* . 2 * 1 > = • 4 
) //y , 

R M -t = ( P M - b « ) / 7 • 

VMC=Nm- H,!/.] 1, 

J=P«-1 }«)/?, 

IF ( fi , |_ J *0. ) ^u:(l, 

IF (^-^.LT.ft. ) = 

IF {^^r.LT.O. )^-ic = o. 

IF I^^O.LT • 0 • ) ^M-) = 0. 

1 F (wma.GT. 1 ♦■) r?M A= j . 

. 1F‘(R*m.oT. 1. >-Mn=l . 

■ IF t^fc.OT.l. >^*C=1 . 

. IF (*Mn.c,T. 1 * > K««|)=1 . 

■*. 1 *<: > =.S.**A»*3 

L^tirO ( « * 4 . d > =* 1 -WMA* . 1 
. ; \ £Nt^< 4*1 • 3)s.<>»>4Mf4». i 

£^h*b s.6-*y ri *.6 . 

t.4F>b (m,4«3) = . I ' . 

ENf "f'-> < :*•» b . 3) = . 1 •£ 

dMtrtb (M* 1 , 4 )= . .-'«■ . 

ENl^=> ( -1.4.4) . b-^MC«* 1 

ENF rf.b ( ■*. * S ♦ 4 ) =. . 1 .»mH* • 0b ♦ 4C* • <?b 
iNE-fb (4*6.4)=. I* xMij* ■ 0 6 ♦ P -*C* . ^b 


tte 
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ENtRia (M* 1 . S> = .:5-RMH*.3 
Enlk^ <M,4*S) = .S-RMtJ«*2S 
£NE*b < m ; 5 . 5 )=*«*» . 5 
ENERb (M,fe,«S) = c?m.-}*.OS 
ENERb < M« 4*6) = , 3-RMH* •3»R«D < > . 1 
ENtR? ( m« 5*6)>, 3 1-kMD*. J5 
ENE^.b <M*6,b) =• j*tfMe*.3-RM0».b5 . 
tNERs ( M ♦ 7 1 6) = . 1 ♦*«[)* , i 
tNERb =RM0».6 
24 0 CON f I N'JE 

C NORMAL lit ENEWb R£ It row E \CH CONSUMER LEVEL 

." DO m=1 /■."•••• 

00 '24 1 L=UO 
i»ecoNjS(M,L)=o. 

DO 241 J=1 ,H 

StC0Ms<M f L>=!>EC0N5(M,L) ♦ ENFRa ( y « J « L ) 

■ 24 1 COnIT INJF 

DO 24? M=1 ,?b 
DO 24? L = 1 .6 
DO 24? J = 1 , M 

ENtR:>(M*J*L>=ENER*(M,j.L>/StC.ONb<M,L> 

24? CONTINUE 

00 21} L=1 «* 

I b < L) = U 

211 CONTINUE • " 

C kFAO IN lM AND F ** FACTOW^,ML,Y AND A GROWTH RAJEi> 

DO 22 1 L = ] .b 

-<EAt)lS,?20) H(L)*I?(L).TJ(L).l4(L).IS(U)*C0NSTlM(L>*FACTt^(L)*ALY 
1 b« < U 

2?0- FORMAT (blb.FlO. 1 tPf 10. S) 

IF (L.GT. 1 ) GO TO 223 
• *RITE <h*2?2> 

222 Fo-vnat ( 4 x ?H I 1 4X2H! ?4X?Hl 34x?ril44X2HlS3x MCONbT IM4X6 hF ACTEmSxSHhl^G 
IV) : ; ' • •• ... - . •■•■■. ■ ' 

221 CONTINUE -•'■-• 

VRI Tt (6.224) II (L) •I2(L)*l3(L)«l4(L>*lb(L> *CONi>TIM <L) .FACTE MIL) »RLY 
1GRIL) 

224 FORMAT (bib, F10. 1 . 2F 1 0 . -> ) 

221 CONTINUE 

DO 226 I=l*?b 
DO 22S L=1 .6 
on ??6 m=i,s 
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Gt^A* ( I * L t M ) s.l * U 
GEPY- < I »L* M) =J .0 
Ft'^AJ ( I ,l.M) =l,(j 
r EPf-c < I ,L*M) =1 .0 

^?S CONT iMOt 

C EXP b^J^TH E ACT O’? FOP YOUNG* l. f > ED* AUJLT F OP SAbELlNt YEaP 
DO <?bb . L = 1 ♦ 6 
»UMtxr,p{L) =0.0 
:>o i = i *?g 

EXPG«t(l*D = <EXP{ .nS« <27-1 ) 1-1. > 

bUMLXOP t L ) = S J*F XGP ( L > ♦►XPoP ( I » L ) 

2?N CONTINUE. 

JO 2 30 1 = 1 *2b 

FfcfN(l«L»l>s tv* r 7«l 1 *L > /PU^tXGP (L ) 

2 30 CONTINUE 

COPP(L) =BLY(w <L> /(EEPY* ( l3,L*l)*l«0> 

COPP ( L ) = (rtLYGP(L)-l . ) /Ffcpr^u 3,L. 1 > 

DO 2bb 1 = W?S 

FEP7P I I »L • 1 ) sFFPYP ( I *L« 1 ) *COpP < L ) *1.0 
FE^Y- U ,L* 1 ) = 1 .0 
"EPA w ( I * L • 1 ) = 1 . 0 
Z*b CONTINUE 
23S CO IT 1 M*Jt 
N =1 
Nl = l 
N2 = b 

2N0 CQ\'T 1 NUE 

C DtTEPMlNc Afi:>Ai_ * A AN I T ON STATION 8Y STATION HASIb 

C WPAPAM-rfATt* PAPAMETE^b-^TAT [ON.MACTIDNG rtQ.IVAP PO ADS , T P I N I T Y PlVFP* 
C HSCOTHE».bALlNI TY (PPT) .TOTAL N I T PObEN ( MO/L ) 

AAoAL ( N) =0. 

A AN I T <N) =0. 

IF (N.GT . 1) GO TO Si M 

C DETEPMINE AAbAL * AANIT ffrOM CHANGE IN F«EbH*ATEP i WASTE OI^CnAPGtb 
DO bi IS 1=1 *?7 

PEAU<S*SUbJ (wPAftAMI J* I ) . J = 1 »7) 

SUb F OPM Ay ( 7 F 1 0 • ? ) 

IF (1 *r,T. 1) GO TO 512? 

WP I Tt { h * 5 l 2b) 

S 1 ?b FOp.MAt ( iX7HbT ATION3X 7HFPACTSP3X7OFPACTTP2X^HFPACTMbCSXSH0THEP7XArt 
IbALSXSHT NIT/) 

5 1 ? 7 CONTINUE 

aPITE (b.Sl 1 b) (wIPAPAM ( J, t ) ♦ J=1 , 7) 
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AA^Al (N) sAAbA|_ IN) *NPARAM (.*>• I )'/2 7« 

AAnJI T(N > = A AM I T < N ) ♦ *P Ak AM ( ( * I ) /d7 • 

5115 CONTINUE 
N Y E A 9 = N 
t>0Sg COMTlWuti 

C *ATE9 \> wA^Tt INI TY ■> I vEP» OTHt-«* TOT Al-wATE 9 CFb,*AbTE 10E06 

C PQJNDb 500^> Y E AU 

-*£40 <S,bl?OM*ATfc9bO<N»M .< = 1«4> , UAbTtbO<N,L) .L=l t4) 

51?0 FOh/MAj <MflO.?> 

• lF<*ATf_PbO(N. 1 ) .EO.O. ) GO TO oil -i 
l F (N . (*T • l ) 00 T O *>128 
rfPlTt (0*5129) 

5129 F09^Aj ( 1 XI ■YHwATEK OI^ChA^GF«C- 0^«lX?9HVKAbTE OlbCrtA^OE MILLION* POUn 
lDb bOO) 

*9lTE(O*5130) 

5130 FOP8ATOX7HT PIvEP 7x3HHbCbx5HOTH£P5xSrtTOTAL3x7HT 9 I v/F 9 7X 3HH5C5X 
15H0THFP5X5HT0TAL ) 

5 1 ?tJ CONTINUE 

30 7003 MM=1,4 

*A:>TE bO (N.NM) =WAbTE:>0 < 1 * N ^ > 

IF(N.FO.l) 00 TO 7003 
IF {.M.GT .5) GO TO 7 0 0 A 
*ATt-<bO (N,N>> = *ATt'9bO ( 1 ,\j*) /i .21 
00 TO 7003 
7004 CONTINUE 

»AT£*bO<N*NM) =1 • l*WATE-<bn (N-l »NM) 

7003 CONTINUE 

rfrtl T£ <b'5]?0) (4ATF.k:>0 ( N*K ) *K=1 «4) • ( MfAbTEbo (N't.) *L>1 «4) 

DO bl T’A < = U4 
9WT U • K J * 1 • 

9<rfb( l,K)=l. 

IF (N.Fu. 1) GO TO 51 39 

9WT(N.K)=^ATEPbO(N,K)/«ATt9b'j(N-l*^) 

' 9Wd(\,k )=wASTtiO(NfK)/*AbTf i>0<N-l,<) 

bl39 CONTINUE 
N = N ♦ 1 

GO TO. 5050 
5119 CONfJMOE 
N = N- 1 
MNN = N 

N=NYLAU g 

51 IB CONTINUE *r 

•iTOhU <N) =0. 



«SQlbU(N)=0« 

DO 5121 1=1.27 

*$Ol < I ) = (wpAPAM<^,I). , *PAPAM(3.I> tt -*rf5(N.l)*rf‘ :> APAN(4.I) 

.rfPAWAM<5. i ) #P*b(N. i) ) 

**T i) i < I ) = ( ^r-A^Avi I ) ♦W^A^AM ( 3. 1 ) *P*T < N * 1 ) .rfPAPAM (4* I ) 

l*PNT(N.2>**PAPA4m,I>* u WT<N.j)) 

rfS;>liO(M=wSD15U(N) ♦ *^01 ( I ) /27. 

*T0 1 P'J(N) = *TD I iO ( N J ♦WTO I ( I ) / df • 

DO *l4l H = 2.5 

IF ( i* • *JE • 2 ) OD TO 5122 

*pArtAvi(M' ] ) = 4.PAi*A4{M' I ) / *.jnl ( I ) 

GO TO 5141 - 

5122 CONTINUE 

WPAHAM (M. I) =*P4PAN (4. 1 ) »U*T (N.4-2) /#TDl (I ) 

SIM CONTINUE 
5121 CONTINUE 

IF (N.eo. 1 ) GO TO SI l 7 

AAnIT(M) sAAMll <N-1 )*«riolSJ <N> / rfTUIbj(N) 

KY/m- *ATtPSO< 1 .4) *AlOG(*a:> 4_ ( 1 >/12,)/AlOG(.5) 

AApal(N)= 32.*(.s-)**(*katf-7sO{\*<*)/XYZW> 

AA5AL<N>=32.* ( AA5AL IN-1> / 32. ) »**IOISU(N) 

**mE<5.S30i > AAaAL ( N) .AASAL^-l) .rfTDlPUlN) , 

1 AANIT (M) » AAM1T (N-l ) t-/^[)lbO 

2<N) *M 

5301 FOPHAt ( bElb.4, I 3) 

5117 CONTINUE • 

SI 25 CONTINUE 

DATA I'M. 1 ) = WATE^30 (N.M 
DATA <N.?> =^A^TE5»0{N.A) 

DA T A ( M, j) =A A 3 AL <N) 

DATA (N,4> =AANIT (N) 
f G TN) = 1 . 1 - • 1 * A AN 1 T < N ) / A An I T ( 1 ) 

PPG(N) =.821 l“.0207 a AANlT (N)— .0129* (AANlT ( N ) ** 2 * ) ♦ • 0 44 9° A Ab AL ( N ) 
NK=\l-2 

lF(NK.Lt.l) NK=1 

DMNG ( N ) = ,094b* (2S.-AASAL INK) 

DATA (N.S) =PPG(N) 

DATA(N,6)= FG(N) 

C CALC OE < 1 J 4 PP<2> 

IF (.M.LF.l ) GO TO 450 

IF (N#F g • 2 ) POGft = PPG(N> /PPblN-1 ) 

IF(N.Gt.?) P»GR = PPG(N-1 > /^G( N-2) 

4s4 continue 
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< = N- 1 

DwF=R*T (N»4) 

DEE = R.Vb(N.4) 

00 453 I=i«?5 

F0U0(M* I 1 1 ) sFOOIXK* I » 1 ) * ( .4?* ( u*F** .50 ) ♦ . 58*PP0R > 

FOOU<N«If.?) =FOOO (K ♦ I * 3 ) * D RG ( N) /°PG < N- 1 ) 

F 000 ( M * I * 3 ) -F 000 (K # I * i ) 

453 CONTINUE 
450 CON T 1 NU t 

00 300 NY=N1.N3 
DO 303 L=1 *6 
TOONS Im<L>=0. 

TCON&EM (L) =0. 

303 CONTINUE 

DO 3b7 1 = 1 *?t> 

00 3b7 L*lt6 
FORE* ( I « L > =0. 

CONbtM { I ,L ) =0, 

CONSIm < I ,L) =0. 

3*7 CONTINUE 

DO 30 3 1 = 1 *36 
00 30 3 L = 1 *6 
IF (NY. or. 1) GO TO 30 7 
RI ( L ) = F L 0 A T ( I- 1 1 (L) ) 

R3 <L)=FLOAT ( I 3 (L ) - 1 1 (D.) 

IF ( R3 ( L ) .LT.O. ) R3 (L) = R3 <L) *3b. 

IF ( R 1 (L) . L T . 0 . ) GO TO 3M 
IF ( r? 1 (l> .EO. 0 . ) GO TO 304 
307 CONTINUE 

i6(UsIb(U*l 

IF(IfXL) .GT.3b> 16 (L) =1 . 

R=I.6<L) 

<N ( I *L ) =N 

C CORRELATE m ANO I FOR EACH L 
GERAR ( I «L* 1 ) =FEPA« (H.Lf i ) 

GEP YH ( I*L. 1 ) =FFRYR(M«L» 1 ) 

C CALC GROwTr, RATE EFFECT OJF TO tNVl* % POLL 
I F ( N.EO.l ) GO TO 388 
C CALC GROWTH rate FOR CURRENT Y f AH 

C CURRENT VEftk GROWTH RATE - F00.> AVAILABILITY AND CONSUMER DENSITY 
FAVRT ( I EL ) =0. 

IH=I-1 

IF< IH.LT'.l ) IH=I 



r> o 


IJ=1d(L)-1 
IFMJ.LT.i > 1 J=1 
FQ < H =400./63b. 

FQ ( ) =400. /5*»* 

FQ ( J ) =400 »/ 1 • 

F(m> =400./750000. 

FQ<5> =400./437S. 

FQ<b) =400./l?425. 

FQ ( n s400./566*R. 

F Q { 8 ) = 400,/3430 7. 

C Oh TERM INh GROWTH CHANGE OUh' TO CONSUMER DENSITY 
GRCOMSO=1. 

IF ( I. E 0.1) GO TO 51 01 
GRCO(1.I.L)=1.0 
GPCO (N*I*L)= • 

1 <W*.0275*I1.-COnML<_> /(C0Nb7L<ltlH*L)*4.n>/ 

? ( 1 .♦.0275<M 1 .-.25) ) 

GPC0N50=GRC0(N« I*L>/GrtCO(N«l*ItL) 

5101 CONTINUE 

00 2r>6 J=lt8 

IF { J.ll.3) GO TO S 1 50 

IF( I .FO. 1) F 000 ( N • I H » J ) = F 0 10 ( N- i ♦ IH, J) 

IF ( N • l.E • 2 ) GO TO 5150 

IF ( 1 .EO. 1 ) FOOD <N« IH, J) =FOO,)(N-1 .?G» J) 

IF II.E:j. 1)F00D(M-1 • lH*.j)=FOOU(N-2,2*fJ> • 

5150 CONTINUE 
FF < J) =0. 

C IF<tNE&MU,J,L) .EG.O.JGO TO 25* 

FF(J)=U.*.01*<1.2*. 75*alOGLO<FOOD(N, IH,J»«FQ(J))))/ 

1 (1. + .U1M1.2*. 75«ALOG10(FOOD<N-1 »IH,J)*F0(J)>)> 

FAV< J)=ENE3SU J, J,L)*FF <J> 

FAvRT ( I .L)=FAV^T ( I tL) ♦FAy(J) 

25* CONTINUE 

WRITE <6*<*70 > FQ < 1 ) ,F()(2) •FQU) ,F0 (4) *FQ(5) ,FQ lb) »FU<7) » FQ<8) *I*L*N 
WRITE (6*470 ) FF (1) ,FF(2) *FF(3),FK(4) ,FF(5> *FF(6)*FF<7) ,FF(8) , I,L,N 
470 FORMAT (HF10. A, 316) 

R=FAVRT ( l.L) 

C TO ALLOW FOR FOOD SUBSTITUTIONS 
IF (R.LT. 1 . >R=R»*.H . 

RR = R 

c growth rates per oay*mooel perjou 14 days 

GEP YR ( I * L * 2 ) - ( RagpconSQ )**!<*• 

F AvR T ( I * L ) =0 • 


10T 




DO 25R J=l,8 
FF ( J I =0 » 

C IF (EnE*S K2h* J»L) .EO.O. > 60 TO 258 

FF(J)s{l.*.01 e U.2-».7B»ALOG10(fOOO(Nf IH,J)*FO(J) )))/ 

,1 <1 *01 * ( i ,2* • 75* AL.D61 0 (FOOD (N-l «IH,J). b F3(J1 J 1 > 

FAV(J)=ENER5(26*J*L)*FF (J) 

F Ay/t^T ( I,L)=FAvRT< IfL> «-FAy( J) 

25B CONTINUE 

R=F A yRT (I ,L) D 

C TO ALLOW FOR FOOD SUBSTITUTION 
IF(*.LT,1.)P=R**«8 

C GROWTH RATES re P DAY,MOOEL PERIOD I 4 DAYS 
GEPAR ( I *L*2) = (R*GRC0NS0)«»1<*. 

C uRITE (6*6O0O)GPCONSn*PR,GE»YWU*-*2) *R*GEPAR( I*Lt2) « I *L 

6000 FORMAT (5E16.H, 21b) 

C oEPAR(I,L*2) = (.3^6*R-.1A21*H**2.-». 0232*R**3 . > /.2555 
GEP Y H(I»L*3) = 1,0 
GEPARf I*L«3) =1 *0 
GERARU*L.3)=FG(N)/FG(N-1 > 

GEPYR U.L*3)=FG<N)/Fu<N-l> 

GEPaR ( l » L , 4 } = 1 . 

GEPYR ( I » L * A J s I . 

C CALC GPOwth RATE EFFECT DUE TD DISTRIBUTION 

C «HEN FPtRH wATEP DECREASES SHRIMP AND CRABS HA yE LESa AREA 
C WHEN FPESh WATER INPUT DECREASES /P PRIM AND MlD CARN HAVE GR AREA 
GEPAR (1,1,4): .500* ,500*AASAl (N) /AASAL (N-l ) 

GE°YR ( 1*1,4) =GEPA»( 1,1,4) 

GEPAP ( i ,3, A) * <!.♦( .54 » < (OMNG(N) /OMNG(N-l) ) -1 • > >> *M 1 ,/2b. ) 
GEPYR( I,3,a)=GEPAR ( I , 3, A ) 

GE°AR( I « A, A) =GEPAR ( I * I *A) 

GEPYP ( I ,A,4) =GEPAR ( I, 1 ,4) 

G£PAR( I ,5,4) sGFPAR ( I » 1 , A ) 

GEPYP ( I ,6* A) =6EPAR< 1,1,4) 

4S2 CONTINUE 
260 CONTINUE 

C STORE PRIOR YEAR NET GROWTH RATE 
FEPAR ( I ,L ,5 )=GEPAR(I,L,S ) 

FEP YR ( I ,L *5 1 =GEPYR ( I *L *5 ) 

?A8 CONTINUE 

. C DETERMINE CORRENT YEAR NET GROWTH RATE 
510 CONTINUE 

IF (N.FO.l) F£RARa,L,5)=GtPAR( I,L» 1 >’ 

IF (N.EO.I) FEPYR( I ,Lt5)=GEPYR(I,Lt i ) 
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GEPY^n ,L.5) =FEPr^( I,L.S) <, r,eprRU,L*2)«GEPr»( I,L,3)*GEPY^(I,L*A) 

GEPAP ( I ,L*5)=FEPAR( I*L*5)*GEPARU,L»2J*GEPAR< I *L • 3 > *GEPAR ( I*L*A) 
(U sPLOAKIfc (L)-l ) 

208 CONTINUE 

RX12=R] (L>/P2(L) 

CONil.M 1 1 *L ) = 0* 0 

IF < P X 1 2 • GT • 1 • ) GO TO 205 

CON*Im( [ ,L ) =*IN (HKD *3,1 4 1G/R2 <D ) *C0nStIN(L> 

205 CONTINUE 

N3(L) =PI <L > -FLOAT ( 13 <U-H <L> ) 

IF(P3(U) .GE.2G. I»KLI =R3(u>-26. 

IF(R3(L).LT.O.) GO TO 210 
(L ) sFLJAT ( I 4 (L ) - I 3 t L > ) 

IF W4(i_) «LE.O. ) R4 (L 1 =R4 (L ) *?b, 

R5(U=Pl(L)-F L 0AT(IS(L)-lHLn 
IF(R5(L) .6E.26. J R5 lL>=Rb (L>-26. 

IF<P3(L) .LT.O.JGO TO 2*0 
R6<U=FL0ATU 1 (L)-l-IS(L) > 

IF (R* (L ) .LE .0. ) R6<L> =R6(L> +2-5. 

RX56=r5 ( L > /R6 (L ) 

IF (PX56.GT . 1 . ) GO TO 2*0 

212 FORl* { I « L ) = YCON5 II (l > * < R5 <L > / < *6 (L ) ♦ 1 . > > 

IF (FORt* ( I *L ) .LT.O. ) FORENt I ,L) =0. 

GO TO 2*1 
280 FORE* U *L> =0. 

2*1 CONTINUE 

PXjH = P3a>/P<Hl> 

IF(RX34»GT*1«) GO TO 2*2 
CONPtM { I ,L ) s 

1 SlN(hf J C u ) ^ 3 . I4 16/P^<L> >*FACTEM lL) ft YC0N5lL (L ) ♦FOR EM ( I *1 > 

GO TO 28 3 

2*2 CONatMUtL)aFOREM(NL) 

2*3 CONTINUE 

GO TO 211 , 

210 CON2>Em ( I * L )■ - 0 * 0 

211 CONTINUE 
A*GPO>V=ACON5IL (L) 

PERGPr (I *L) = YCONGIL (U* (GEpyr( I,l*5)-K) 

IF (PfcRGWY (,I*L) • L T » 0 » ) PEPORY ( I *L > =0 . 

3 ERuR A ( { » L) =ACONG II <L ) * (GE^AR ( I *^«5)-l « ) 

IF ( * • EG • 1 )CON5tMCI,L}=YCON5U{U«-GEPYR(I,Lt5) 

ACON^IL(L) =AC0N5IL (L)*GEPAR( I,L«5) /F£PAR( I,L*5) 

IF(rt.FO.i) ACON^EM(U =ACONi>0*(LJ 



noon 


RESULT ( N *L 1 3 ) = ACON^EM (L ) 

YCONblL <L) =YCONblL<L> *GEPYR ( 1 *L*5) ♦ CONS I M { I ,L> -CONSEM ( I ,[_> 
IF<H.P(3.1 ) ACONbOM< t)-U. 

C ADULT EH=GR0*TH-CHANGE IN BIOMASS IN SYSTEM 

PERAE M 1 1 *L > = RE«GRAU*L)-( AC0NSIL (L)-ABGROW) 

ACONbOM (L)=ACONSOM(L) ♦ PE°AEM(I,l) 

CONb£H{ I*L) =CONb£M( I,U ♦PEPAtMdtL) 

GO TO 202 

204 ACONbJL (L) =S ( I ,L ) 

YCONB II ( L > = 0 • 

CONblM ( I ,L ) =0 , 

CONbEM(I,L)=0. 

16 (L ) = 1 

202 CONTINUE 

CONblL < L > ^ACONblL (L) ♦YCONb I L ( *_) 

C0NS1L(I«U =AC0NSIL<L) 

IF ( YCONb IUL ) .LE.l.) YCONbIL (.) = 1 . 

C0NS2L U *L ) sYCONS IL (L ) 

C0N>3U( ItU=CONSlL(L) 

C STORE CONSUMER BIOMASS LEVELS BY YEAR 
CONS7l<N*I*L>=CONSIL (L) 

RCALbAM(I,L)=CONS3L { I*L)/S <I,L) 

RCALbAM BASED ON YOUNG BIOMASS 
XYZ=b ( I,L>-C0NS1 l{ I,L) 

IFIXYZ.LT. l.)XYZ=l. 

RCALSam(I.L) =CONS?l ( I.D/XY2 
IF-IN.GT.l >PCALSAMU»U=C0NS3 l(I*l)/CONS 6L( 1,0 
IF (L.GT.b) GO TO 254 

ll=lo 

FOOLMN, I.LL)=CONS3L (I*U 
2S4 CONTINUE 

TCONS|m<L)=TCONSIH(L) ♦COnSim ( I ,L> 

TCONbFM{L>=TCONSEM(L) ♦CO nSEM U,L) 

BEbULT (N.L. 1 ) =TCONblM (L ) 

RESULT <N.L *2) =TCONSEM(L) 

REMIM=tCONSEM ( l ) /TCON blM (L) 

RESULT (N*L ,4) =REMIM 

203 CONTINUE 
IF <NY„LT.N2) GO TO 200 
DO 2Bft 1=1,26 

np=np ♦ 1 

XP (MR) sFLOAT (NP ) /26. 

00 6100 L- 1.6 . 
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rcow^l (NP*L)=C0NS3L ( I.L) 

6100 CONTINUE: 

IFU.Gt.I) GO TO 266 
-(RITE (6,269) 

269 FORMAT < 6 X 5 HF 0 OrH 5 xSMFO 0 D? 5 xSHF 00 Q 35 x 5 HE 00 D 45 X 5 HF 00 DSSXSHF 0 DO 65 x 5 HE 

lOOO/^XSHFOOOBGXlHI ) 

268 rfRUE ( 6*270 ) FOOD <N* I. 1) * FOOD ( N* I *? > *F00D ( N . I . 3 > * FOOD f N , I , 4 ) , FOOD <N 
1 , 1 * 5 ) » FOOIMN, 1,6) , FOOD ( N , 1 , 7 ) , FOOD ( N, I , 8 > , I 

270 FORMATUx*BF10.1,I6) 

298 CONTINUE 

DO 6110 L= 1 , 6 

C SKIP DETAILED PRlNTOUTb-SlOO 
GO TO Si 00 
*RITt (6*219) 

219 FORHAj (9X7HTC0NSlM3x7MTC0NSEM‘*X2HNYSXlML5x5HREMlM3x7rlAC0NSEM) 

-(RITE ( 6 * 2 lRJ TC0NSIH<L » *TCOnSEM(U *NY,L*REMlN,ACONi>EM(U 
218 FORMAT <1X,?F10. 1*216, Flo* c**F10«1/) 

DO 310 1=1 ,26 . 

IF(I.GT.I) GO TO. 311 
-(RITE<6,300) 

300 FORMAT (Sx6HGEPYm4X6H&EPrR2Ax6nGEPYR34X6HGEPYE44X6nGEPYR55xlrtL) 

311 CONTINUE 

RR1TE (6*30 1 > GEPYR ( I ,L, 1 ) , GEPYR ( I ,L »?) * 6E P YR ( I * L * 3 > , GEPYR ( I , L , 4 > , 

1 GEPYR ( I *L *5) , L 

301 FORMAT ( 1X.SF10.4, 16) ' '• 

310 CONTINUE 

DO 312 1 = 1 *26 

IF ( I .GT. 1 ) GO TO 313 

-(RITE (6,302) 

302 FORMAT (5x6HGEPAR14x6HGEPaR?4x6H&EPAR34x6HGEPAE44x6HGERARSSX1HL) 

313 CONTINUE 

-(RITE (6*301 ) GEPAR ( I ,L. 1 ) *GEPAR ( 1*1*2) , GEPAR ( I ♦ L * 3 ) *GE p AR ( I *L *9) , 
1GEPAk(I,L*S)*L 

312 CONTINUE 

DO 216 1=1,26 
IF(I.GT.l) GO TO 215 
-(RITE (6,214) 

214 FORMAT <aX7HAC0NSIl3X7HYC0NS!L3X7HTC0NSIL3X7HRCY/PRY4X6HCDNSIM4X6HC 
13N3EM4X2HNY5X1HI5X1HL5X1HM5XSHF0REM4X6HGEPYRSAX6HGEPAR5) 

215 CONTINUE 

-(RITE (6,21 7 ) CONS 1 L ( ItL > ,COnS2l ( I ,L) ,C0nS 3L ( 1*0 *RCAl_SAM (1,0 »CONSl 
1 M ( I , l ) *C0NSEM(I,U ,NY, I,L *KN(I«U) *EOREM( I,|_> , GEPYR ( I,l,5) , GEPAR (I 
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2*L*5) 

217 FORMAT UX*6F) 0 . 2 • 4 1 6 * F 1 0 .2 *2F 1 G . 4) 

216 CONTINUE 
5100 CONTINUE 
5110 CONTINUE 

C DETERMINE FOOD COEMPTION by TYPEtCONV EFF=.l 
DO 500 J*1,M 
TAFCDn5(J)=0. 

AVGF0nO( J)sO. 

DO 500 1=1*26 
TPFCONSU,J)=0. 

500 CONTINUE 

DO 504 J^l .8 
DO ^>01 1 = 1*26 
DO sU3 L= 1 *6 

^FCONb< I, J*L) =(PERORY ( NL)«ENtRb(I.J,L)* J> ERGRA(IfL)*ENER^I26* J.L> > 
l/.l 

TPFCONbd . j)=TPFCON^( I ♦ J) ♦°FCDN6 { I * J*L > 

503 CONTINUE 

TAFCON^I J) =TAFCOnS(J) ♦TPFCON^C I* J l 
AVGFOODI J)=AvGFOOU( J> ♦FODD(N*I *J)/26. 

501 CO NTTNUE 

RFCA ( J) =TAFCDnS ( J) /AVGFOOD( J) 

504 CONTINUE 
rfRITE (6* 505) 

505 FORMAT (3XRHRFCO/ AVI 2X«HFFCO/An22X8HRFCO/Av32XBHFFCO/A\/42X8HRFCO/Av 
152XdNRFCO/Av62X8HRFCO/AV72XBHRFCO/Av85XirtN) 

rfRlTE (6*506) RFCA( I ) ,RFCA(?) , RFC A ( 3), RFC A (4 > , RF C A ( b ) * RF CA ( 6 ) *RFCA <7 
I > *RFCA ( 8 ) «N 

506 FORMAT <1X*8E10«2*I6> 

200 CONTINUE 

c store pr year food and biomass levels - * 

DO 2b3 1=1*26 

00 265 L = 1 * 6 

C0N^4|_ ( 1 .1 ) =C0NS1L ( I *L) 

C0nS5L{I*L)sC0nS2l(I*U 
CONSbL ( I *L ) =C0NS3L ( 1 *L ) 

265 CONTINUE 
263 CONTINUE 

C SUM TOTAL BIOMASS EM FOR YEAR 
RESULT < N * 6 * 5 ) =0* 

DO 490 L=2* 6 

RESULT (N,6,5) =RESULT (N,6, 5) *REPULT (N,L*2) 
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490 CONTINUE 
M*N'+ i 

N1=N 

N2=N1 

IFlN.GT,NNN) GO TO 4?6 
GO TO ?90 
476 CONTINUE 
N=N-1 

00 480 L = l*6 
DO 484 K = } «N 
IF (K.GT. I) GO TO 4R2 

481 FO^rtAy ( 7X4HAF* r V7X3HA*D5X5rtAA:>AL5X5HAAN IT7X3HTI97X3HT£M3 x7hAC0NSEM5 
1x5UREmIN5x1HNSx1Hl7x3HPPGBx2mFG5xSHTEmYR) 

482 CONTINUE _ , , . 

WRITE (b* 48 3) DATA (Ki 1 i * OAj A <K ♦ 2 ) »OAT A (K« 3 ) * DATA (K *4) * RESULT CK*L*JI * 

1 RESULT (K .L *2 J « RESULT (K «L « 3 ) t RESULT (K*L*4)»K*L *DAT A<K*5)*DATAlK»fe) 
2*R£sULT(<«6,5> 

483 FORMAT (1X»4E 10. 2* 3£ 10.2* E 10*2*216* 3EIQ.2) 

4K4 CONTINUE 

DO 465 NX s 1 * NP 
TCONS(mx) =YC0NSI (NX*L> 

485 CONTINUE 

IF(L.EO.l) YM=13H2 OOPLANKt 0N 
IF (L.E0.2) YM = 13HHER8lV0RE 
IF (L.E0.3 ) Ym=I3hOMNIVORE 

IF(L.F0.4) YM=13hPRIM CAWNlV 
IF (L*F0.5>YM~1 3r»M I D CARNIV 
IF(L.EQ.b)YM=13HTOP carniv 

CALL DOIPLT (1 • IN»NP«XP, YCONS , 0 . ♦ 0 • f 0 » * 0 . ♦ 1*XM« 2*YM,11) 

480 CONTINUE 

XM= 1 UHYEAR 
Y8=I3HHI0MASS 
DO 7001 <L=l«9 
DO 7000 K=UN 
ZP(K)sFLOAT <K) 

IF(KL.EO-l) YRESU {K ) =UAT A ( X * 3 ) / 1 0 * 

IF (KL.EO .2) YRFbU (K)=OATA (K»4> 

IF(KL«£Q*3) YRF.bU(K)= RESULT (K ♦ i ♦ 2 ) /RESULT < 1« 1 * 2 ) 

IF < KL . EO • 4 ) YRFSU(«)= RESULT <K«2*2> /RESULT < 1 *2*2) 

IF (KL.EU.5) YRESU(K) = RESULT (K<3*2) /RESJLT ( I *3*2) 

JF(KL.FQ.b) YRESU(K>= RESULT (K>4»2)/REsULT U *4,2) 

IF (KL •E'U • 7J YRESU(K)= RESULT (K *5*2 ) /RESULT 1 1 * 5*2). 


o 
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1F<KL.EQ.8) YRE=>U<K> = RESULT (K.6. 2) /RESULT (1 *6.2) 

IF(KL.E0.9) YRESU(K)= RESULT CK*6.S)/RESULT (1 .6,5) 

7000 CONTINUE 

IEC = 0 

IFtKL.EU.9) IEC=1 

ISVH=kl 

CALL no I PL T < IEC. IN.N.ZP.YRFSU,Q.,0.«G.*2.0.1 .XM.2.YM.ISYM) 

7001 CONTINUE ' 

XM=10hPERIOD 
YNsIOHFRACT food 

DO 7012 -LSI *6 
JO 7013 J= 1.8 
DO /0 1 4 4=1,26 
TRE^U(M) =ENERS(M. j»l> 

ZP{M)=FLOAT(M) 

7014 CONTINUE 
4=26 
1EC = 0 

IF ( J.EO.b) IECsl 
ISYR-J . 

CALL 00 1 PL.T ( I EC. IN.N.ZP. YRFSU, 0 ,.0»*0,*1 •0.1.XM.2.Y4, I5Y4) 
7013 CONTINUE 
7012 CONTINUE 

YM= I 3HBOO. LB/OA Y 
DO 7020 4=1 ,26 
YR£3jm)=FOOD( l .M, 1 ) 

7020 CONTINUE 

4=26 

IEC= 1 
ISYM=1 

CALL DilIPLT <IEC. IN,4.ZP,Y*E*U,0.*0. »0..0. . 1 , XH * 2* YM. I SyH.) 

YM=13HBI0MASS 

DO 7021 L = 1.6 

DO 7022 4=1,26 

YRE^U(M)=AL0G10(b(M,L) ) 

7022 CONTINUE 
4=26 

IEC = 0 

IF(L.Eu.b) I E C = 1 
I sy H=L 
A I =6 • 

CALL OnlPLT(IEC. lN.M»ZP/YREiU.0..0..0..Al, 1.X4.2.YM.ISYM) 

7021 CONTINUE 



’roiJClBILITY OF IH0 

!V4P.T? TC? 


IEC = 1 
I Sy*= l 

00 7023 H=l,26 
YPESO (M > =FOOO 1 1 • M * 7 > 

7023 CONTINUE 

CALL DOlPLT ( IEC* IN*MrZP*Y*ESU*Q.,0.«0.,0. « 1 Y4« ISYM1 

STOP 

END 


11-5 



source water 

FRACTION* 

FOR GALV 

bay stations (copeland 

AND FRUH« 

1970) 

station 

fractbr 

FRACTTR 

ERACTMiC 

OTHER 

SAL 

T NIT 

1 . 

.74 

.07 

.14 

.05 

27. 

0. 

2. 

.50 

.0 

. .10 

.40 

25.22 

• 75 

3. 

.50 

.0 

.50 

.0 

22.8 

• 60 

4. 

.40 

.10 

.50 

.0 

17. 

; l . 

5. 

.25 

• 10 

.65 

.0 

15.5 

1.2 

12. 

.70 

.0 

.0 

.30 

22.5 

.8 

13. 

.30 

. .0 

.0 

.20 

29. 

• 4 

14. 

.75 

• 0 

.0 

.25 

24. 

• 6 

15. 

.47 

.16 

.33 

.04 

19. 

.92 

16. 

.49 

.09 

.18 

.24 

21.3 

.82 

17. 

.50 

.09 

.18 

.23 

19.5 

1.05 

13. 

.43 

.16 

.37 

.04 

'13.5 

.89 

1*. 

.23 

.16 

.55 

.01 

11.5 

1.15 

20. 

.24 

.14 

.61 

*05 

11.5 

1.70 

21. 

.13 

.09 

.78 


14.0 

3.0 

22. 

.26 

.15 

.58 

.01 

14.8 

1.81 

23. 

.2<* 

.42 

.36 

.0 

12.8 

1.90 

24. 

.20 

.54 

.28 

.0 

4.5 

1.00 

25. 

.16 . 

.63 

.22 

.0 

2.5 

1.00 

26. 

.21 

.52 

.29 

.0 

7.4 

.70 

27. 

.23 

.48 

.31 

.0 

7. 

.9 

23. 

.31 

.34 

.37 

.0 

9. 

• 95 

29. 

.37 

.27 

.37 

.0 

11. 

.90 

30. 

.35 

.27 

.38 

.0 

11. 

.8 

31. 

.51. 

.09 

.17 

.23 

22.5 

.70 

32. 

.49 

.09 

.17 

.25 

21. 

.70 

36. 

.14 

.0 

• . 36 

.0 

13. 

5. 
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